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Nonequilibrium statistical mechanics of weakly stochastically 
perturbed system of oscillators 

A. Dymov*t 


Abstract 

We consider a finite region of a d-dimensional lattice, d £ N, of weakly coupled harmonic oscil¬ 
lators. The coupling is provided by a nearest-neighbour potential (harmonic or not) of size e. Each 
oscillator weakly interacts by force of order e with its own stochastic Langevin thermostat of arbi¬ 
trary positive temperature. We investigate limiting as e —>- 0 behaviour of solutions of the system 
and of the local energy of oscillators on long-time intervals of order e^ 1 and in a stationary regime. 
We show that it is governed by an effective equation which is a dissipative SDE with nondegenerate 
diffusion. Next we assume that the interaction potential is of size eA, where A is another small pa¬ 
rameter, independent from e. Solutions corresponding to this scaling describe small low temperature 
oscillations. We prove that in a stationary regime, under the limit e —► 0, the main order in A of 
the averaged Hamiltonian energy flow is proportional to the gradient of temperature. We show that 
the coefficient of proportionality, which we call the conductivity, admits a representation through 
stationary space-time correlations of the energy flow. Most of the results and convergences we obtain 
are uniform with respect to the number of oscillators in the system. 
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1 Introduction 

1.1 Motivation and set up 

Investigation of the energy transport in crystals is one of the central problems in nonequi¬ 
librium statistical mechanics. In particular, a derivation of the Fourier law and the 
Green-Kubo formula from the microscopic dynamics of particles is of great interest (see 
[Pei, BoLeR]). In the classical setting one investigates the energy transport in a Hamil¬ 
tonian system, interacting with thermal baths of different temperatures. This interac¬ 
tion is weak in geometrical sense: the thermal baths are coupled with the Hamiltonian 
system only through a boundary of the latter. Unfortunately, for the moment of writ¬ 
ing this problem seems to be out of reach, due to the weakness of the interaction. In 
this setting even the existence of a stationary state in the system is not obvious (see 
[EPRB, EH, RBT, Car, HM, CEPo, CE]) 1 . That is why usually one modifies the sys¬ 
tem in order to get additional ergodic properties. Standard ways to achieve that are 
i) to consider a weak perturbation of the hyperbolic system of independent particles 
([DL],[Ru]); ii) to perturb each mode of the Hamiltonian system by sufficiently strong 
noise ([BoLeLu, BeO05, BaBeO, BoLeLuO, LO, BeKLeLu, BeOll]). Following the sec¬ 
ond way, usually one perturbs the Hamiltonian dynamics by very special noise, preserv¬ 
ing the energy of the system ([BeO05, BaBeO, LO, BeKLeLu, BeOll]), or one couples 
each particle of the Hamiltonian system with its own stochastic Langevin-type thermo¬ 
stat ([BoLeLu, BoLeLuO, BeKLeLu]). In the latter case the stochastic perturbation has 
rather natural structure, however the energy of the system is not conserved. 

Let us discuss some of results mentioned above in more details. In [BeO05] the authors 
study a chain of harmonic oscillators, perturbed by energy preserving stochastic exchange 
of momentum between neighbouring oscillators. Coupling the first and the last oscillator 
to thermal baths of different temperatures 7i and 72, they prove the Fourier law, i.e. that 
the stationary averaged flow of energy is approximately proportional to (71 — T^/N for a 
large number of particles N and a small difference 71 — 72• In [BaBeO] the authors con¬ 
sider a system of oscillators, where each oscillator is perturbed by a stochastic dynamics, 
conserving momentum and energy. They study existence of the Green-Kubo conductivity. 
In [BoLeLu] the authors consider a system of harmonic oscillators, where each oscillator 
is coupled with a stochastic Langevin thermostat. Fixing temperatures of the "exterior" 
left and right thermostats they find and study the self-consistent temperature profile for 

1 See also [Tr] and [Dyml2] for a similar problem in deterministic setting. 
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the "interior" thermostats, i.e. a temperature profile such that in a stationary regime 
there is no average flow of energy between the "interior" thermostats and the system of 
oscillators. Then they prove that with respect to this temperature profile the Fourier law 
and the Green-Kubo formula hold. 

In all the works listed above the stochastic perturbation is of order one. It is natural 
to study the case when it goes to zero. Such situation was investigated in [BaOS, BeHu, 
BeHuLeLO]. In [BaOS] the authors study the FPU-chain, where the nonlinearity is re¬ 
placed by an energy preserving stochastic exchange of momentum between neighbouring 
nodes. They investigate the energy transport under the limit when the strength of this 
exchange tends to zero. In [BeHu] the authors consider a chain of harmonic oscillators, 
where each oscillator is weakly perturbed by energy preserving noise of size £ and anhar- 
monic potential of size u, u < e. They investigate limiting as £ —?• 0 behaviour of an upper 
bound for the Green-Kubo conductivity. In [BeHuLeLO] the authors study an infinite 
chain of cells weakly coupled by a potential of size u, where each cell is weakly perturbed 
by energy preserving noise of size e. Expanding the Green-Kubo conductivity in a formal 

series k(z/, e) — i/ n K n (e), they find the term and investigate its limiting behaviour 
n> 2 

as £ —?• 0. It is argued that under the formal limit v —)■ 0 the Fourier law holds with the 
conductivity ^(e). 

In the present work we consider a system of oscillators, situated in nodes of a bounded 
set C C d G N, and interacting via a nearest-neighbour type potential. The system is 
given by the Hamiltonian 

2 

fC(P, 9) = £ (t + u Ai,))+t a £ V( qj ,q t ), (1.1) 

jeC j,keC:\j-k\=l 

where (p, q) = ( Pj,qj)j & c £ v > 0 and for j G 7L d by \j\ we denote the /i-norm. 

The interaction potential V satisfies V(x,y) = V(y,x), is C 2 -smooth and has at most 
a polynomial growth at infinity (for precise assumptions see Section 1.4). The pinning 
potentials Uj will be specified later. Following the second strategy above, we couple each 
oscillator with its own stochastic Langevin thermostat of arbitrary temperature Tj by 
coupling of size £, where the temperatures satisfy 

0 < Tj < C < oo, (1.2) 

with a constant C independent from j. The resulting system obeys the equation 

j t ( h = Pj, j f Pj = - d qj u Mi)~ v d v V ^vqk)-ep j + V2£T j j t Pj, j e C, (1.3) 

k:\k-j\=l 

where /3 = (/3j)j & c £ is a standard |C[-dimensional Brownian motion. Our goal is to 
investigate limiting as £ —> 0 behaviour of solutions and of energy transport for eq. (1.3) 
on long time intervals and in a stationary regime. So, as in the classical setting given 
above, we are interested in the case of weak interaction between the Hamiltonian system 
and the thermal baths, however the weakness is understood in a different, non-geometrical 
sense. Unfortunately, we will have to assume that the coupling constant v goes to zero 
simultaneously with e, with some precise scaling which will be discussed later. 

Our system admits a clear physical interpretation, so not only it can be considered as 
a toy model for a purely Hamiltonian situation, but it is of independent interest. Indeed, 


3 



since the stochastic perturbation is provided by the weak coupling with the Langevin 
thermostats, one can think about the transport of energy through a crystal plugged in 
some medium of a given temperature distribution and weakly interacting with it. However, 
since the energy of the system is not conserved, in our case non-Hamiltonian effects seem 
to be stronger than in the systems above, perturbed by vanishing energy conserving noises. 
In order to have some kind of energy preservation, one could consider the self-consistent 
temperature profile, if the latter exists. 

Since the systems of statistical physics are of very large dimension, it is crucial to 
control the dependence of the systems on their size. Let us emphasize that most of 
results and convergences we obtain are uniform with respect to the choice of the set C (in 
particular, with respect to its size), so they satisfy this physical requirement. 

We started our study in [Dyml4], There we considered a system of rotators, where the 
uncoupled rotators were assumed nonlinear, while the coupling constant v was chosen as 
v = e a , a > 1/2. In the case 1/2 < a < 1 the rotators were required to have alternated 
spins (in such a way we excluded from the system too strong resonances). Using the av¬ 
eraging method of Khas’minski-Freidlin-Wentzell type in the form developed in [KP], we 
found an autonomous (stochastic) equation for the local energy which described its lim¬ 
iting (as £ —> 0) behaviour on long time intervals of order £ _1 and in a stationary regime, 
uniformly with respect to the choice of the set C. However, it turned out that this equa¬ 
tion did not feel the Hamiltonian interaction of rotators, so that in the considered system 
there was no Hamiltonian energy flow under the limit e —y 0, neither on time intervals 
of order £~ 1 nor in the stationary regime. The reason is that under natural assumptions 
resonances of the system of uncoupled nonlinear rotators have Lebesgue measure zero. A 
weak coupling does not change the situation significantly, so that the noise pushes out 
rapidly the perturbed dynamics to nonresonant area, where the Hamiltonian terms aver¬ 
age out. In [FW06] 2 similar results were obtained for a system similar to (1.3), where the 
uncoupled oscillators were nonlinear and v — £ (however, no information was provided 
about dependence of the results on the choice of the set C and behaviour of the system 
in a stationary regime). 

In order to observe the Hamiltonian energy flow, in the present paper we study the 
situation when the uncoupled oscillators are linear and have the same frequencies, so that 

q 2 

u Mj) = y for a]1 3 e C ■ ( L4 ) 

In this case eq. (1.3) is strongly resonant, see below. We choose the coupling constant v 
as 

v = Xe, where 0 < A < 1, (1.5) 

since we are not able to treat the case when £ = o(i'), while the choice (1.5) already 
provides nontrivial results. In the first part of the paper, devoted to the study of the 
limiting as £ —» 0 behaviour of solutions of eq. (1.3)-(1.5), we assume that A is a fixed 
number satisfying 0 < A < 1. In the second part, treating the limiting (as £ —y 0) 
behaviour of the energy flow, we assume that A is a small parameter, independent from 
£ (but usually also fixed). Eq. (1.3)-(1.5) with A 1 describes small amplitude low 
temperature oscillations, see Appendix B for explanation. 

2 See also Chapter 9.3 in [FW12]. 
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Since A will be mostly fixed while £ will go to zero, natural time scale for eq. ( 1.3)-(1.5) 
is f ~ e~ x . It is convenient to introduce a slow time r = et, so that r ~ 1. Then eq. 
(1.3)-(1.5) takes the form 

Qj = e-'Pj, Pj = - A ^2 d qj V(qj,q k ) - p 3 + y/Wjfc, j e C, (1.6) 

k:\k-j\=l 

where the dot denotes the derivative with respect to r. It is well known that eq. (1.6) has 
a unique solution, defined for all r > 0 (see [Khal2]), and is mixing 3 (see Appendix A). 

1.2 Main results: the limiting dynamics 

Our first goal is to investigate the limiting as £ —» 0 dynamics of eq. (1.6) for fixed A 
satisfying 0 < A < 1. Introduce the action-angle variables for the system of uncoupled 
oscillators, given by the Hamiltonian H°, 

2 2 

I = ( Ij) je c e ® |c| , Ij '■= Pj 2 and p = {(Pj)jec e T |C| , p 3 := arg (pj + iq 3 ), (1.7) 

where % denotes the imaginary unit and we put argO := 0. Then eq. (1.6) takes the form 

ij = • • •, & = + • • •, ieC, ( 1 . 8 ) 

where the dots stand for terms of order one (as £ —y 0). Eq. (1.8) describes a strongly 
resonant fast-slow system: the actions change slowly while the dynamics of angles can be 
decomposed to a fast rotation along the diagonal 1 = ( 1 ,..., 1 ) of the torus , and a slow 
drift in \C\ — 1 ortogonal directions £ m , • 1 = 0, 1 < m < \C\ — 1 . A usual way of studying 

of the limiting as £ —» 0 dynamics of such systems is the method of resonant averaging 
(see [AKN]). The latter suggests to pass to the variables (A, -0, 77 ) G x T c l, where 
( 0 , 7 ?) = (- 01 ... , 0 |c|_i,? 7 ) G are coordinates of the vector of angles p with respect 
to the basis (£ x ,..., £^ _ 1 , 1 ), and Tl c l denotes the torus Tl c l expanded in a suitable way. 
Then eq. (1.8) takes the form 

ij = ..., Tpm = ..., p = £ —1 + ..., where 1 < m < \C\ — 1. (1.9) 

Eq. (1.9) is a fast-slow system with the unique fast variable 77 , so the limiting as £ —» 0 
behaviour of the slow variables (/, 0 ) is expected to be given by the averaging of (/, 0 )- 
equations with respect to 77 , which further on we call the resonant averaging. However, 
in order to establish this we can not use the usual scheme which suppose to work in the 
action-angle variables, since writing the yj-equations in more details we find p 3 = £~' +/~ 1 • 
..., j G C, so that they have singularities at the locus {/ = ( Ij)jeC '■ 0=0 for some j G 
C}. Moreover, the resonant averaging of the /-equations has weak singularities there. To 
overcome this difficulty we use a method developed in [KuklO, Kukl3] for a non-resonant 
case and in [KM] for a resonant one. It suggests to find and study an effective equation, 
written in the regular (p, g)-variables, which governs the limiting as £ —> 0 dynamics of 
the (/, 0)-components of solutions of eq. (1.6). To this end, in a general situation one 
should use that the (/, 0 )-components of solutions of the effective equation have to satisfy 
the resonant averaging of the (/, 0)-equations from (1.9). However, following [KM], we 

3 I.e. it has a unique stationary measure and its solutions weakly converge as r ^ oo to this measure in distribution. 
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note that in our case, where the system of uncoupled oscillators is linear, there exists a 
simple way to find the effective equation without using the (/, ^-coordinates, so that the 
singularities does not appear at all. It is based on a passage to a fast-rotating coordinate 
system, where eq. (1.6) has uniformly in e bounded coefficients, see Section 2.2. The 
effective equation has the form 

(Z ; . = Aa Pj ir"-| + /p‘, Pj = -H,H r “-j+VTiPh jec, (1.10) 

where (/3 j)k=i, 2 ,jec are standard independent Brownian motions, while H res (p,q ) denotes 
the resonant averaging of the interaction potential 1 V(qj,qk) (see (2.35)). 4 5 

\k-j\ = l 

Eq. (1.10) has a unique solution, this solution is defined globally, and (1.10) is mixing 
(see [Khal2, Ver87, Ver97]). Let us fix 7 > 0 and some "not very bad" random initial 
conditions (p 0l q 0 ). 5 We prove 

Theorem 1.1. Let (p £ ,q £ )(r) and (p,q)(r) be unique solutions of equation (1.6) and the 
effective equation (1.10), satisfying V(p £ ,q £ ){ 0) = V(p,q)( 0) = T>(po, Qo)- Then 

T>(l(p £ ,q £ )(-)) —^ T>[l{jp, q)(-)) as £ —» 0 on C([0, T], M^) uniformly in C. (1.11) 

Theorem 1.1 is a corollary of Theorem 2.4. This is explained in Section 2.3. Here 
T>(-) denotes the distribution of a random variable. By saying that a convergence is 
uniform in C, we mean that it holds uniformly with respect to the choice of the finite set 
C C 7L d . In particular, with respect to the number of nodes \C\. Uniformity in C of the 
weak convergences of measures through all the text is understood in the sense of finite¬ 
dimensional projections. It means that for any finite set of nodes A C Z d , projections of 
considered measures to A converge with a rate, independent from the choice of the set 
C 3 A. For instance, for convergence (1.11) this means 

E f(l(p £ , q £ )(•))—>• E/(/(p, g)(-)) as £ —» 0 uniformly in C, 

for any continuous bounded functional / : C([0, T], ) H> M, satisfying Supp/ C A, 

where Supp / is defined in Section 1.4, Agreements.6. Here and further on E denotes the 
expectation. 

Remark 1.2. In a similar, but technically more complicated way, it can be shown that 
the limiting (as £ —)■ D) behaviour of joint actions I and slow components on angles if is 
also governed by the effective equation, see [KM], Since it is not needed for study of the 
limiting transport of energy, we skip it. 

Let p £ be the unique stationary measure of eq. (1.6). Each limiting point (as e —> 0) 
of the family of measures {/A, 0 < e < 1} is an invariant measure of eq. (1.3)-(1.5) with 
£ = 0. This equation describes a system of uncoupled harmonic oscillators, so it has 
plenty of invariant measures. The next theorem ensures that, in fact, only one of them is 
the limiting point and distinguishes it. 

Theorem 1.3. Let p £ and p be the unique stationary measures of equation (1.6) and the 
effective equation (1.10) correspondingly. Then p £ p as £ —)■ 0. Linder some additional 
assumption this convergence holds uniformly in C. 

4 Despite that the interaction potential does not depend on the moments its resonant averaging may depend on them. 

5 By "not very bad" we mean that the components have finite exponential moments which are bounded uniformly 

with respect to j £ C and with respect to the choice of the set C C See ass. HI in Section 1.4. 
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See for details Theorem 2.8. Roughly speaking, the additional assumption above de¬ 
mands the effective equation for an infinite system of oscillators to have a unique station¬ 
ary measure, in some class of measures. In particular, it holds if A is sufficiently small 
and second partial derivatives of the interaction potential V are bounded. 

To establish the uniformity in C of convergences from Theorems 1.1 and 1.3 we employ 
a method developed in [Dyml4], Its main tool is Lemma 2.1, where using Z 2 -norms with 
exponential decaying weight, we obtain uniform in C, e and r estimates for solutions of 
equation (1.6). 

In Appendix A we prove that rate of mixing of equation (1.6) is independent from e. 
Jointly with Theorems 1.1 and 1.3 this implies the following result. 

Theorem 1.4. (i) The convergence T> (/ (p £ ,q £ )(r )) -©(/&>,?)( t)) as £ —> 0 provided 
by Theorem 1.1 holds uniformly in r > 0. 

(ii) We have 

lim lim V(l(p £ ,q £ )(T )) = lim limX?(/(p e , g e )(r)) = lim T>(l(p £ , g e )(r)) = II/*/i, 

£—^0 T—KX) V ' T—>OO £->0 V ' ' £—^0 v ' ' 

T—>■ OO 

where p is the unique stationary measure of the effective equation (1.10) and II/ denotes 
the projection to the space of actions. 

For details see Theorem 2.9. Define the local energy of a j-tli oscillator as 

£j(p,q) : = ~ ^ ^ V i ( lT ( lk), so that H v = s ^£ v j . (1.12) 

k£C:\j-k\=l jec 

Let ( p £ ,q £ )(r ) be a solution of eq. (1.6). Since £f = Ij + 0[y) 1 we obtain the following 
result whose full version is given in Proposition 2.10. 

Proposition 1.5. The limiting (as £ —> behaviour of the vector of local energy 

(£j{p £ ,q £ ))jec coincides with that of the vector of actions I(p £ ,q £ ), while the latter is 
described by Theorems 1.1-1.4- 

1.3 Main results: the energy transport 

Our next goal is to study the limiting (as £ —> 0) dynamics of a stationary Hamiltonian 
energy flow in system (1.6). We will assume that A is sufficiently small. 6 

The local energy £f changes due to the Hamiltonian coupling of the j-th oscillator 
with neighbouring oscillators and with the j-th thermostat. Applying the Ito formula, we 
get 

£f — £— * {£), H v } + (thermost. term)^ = ^ ^ J k] + (thermost. term) ? ., (1.13) 

k:\k-j\ = l 

where {•, •} denotes the Poisson bracket, (thermost. term)^ := —p 2 - +Tj + ^2Tj$j, and 

Jkj{p,q ) := { — + 2 Q] - Pk ^ qk :V(qj,qk )} = Pkd qk V(q k ,qj) — Pjd qj V(q k , qf) (1.14) 

is the Hamiltonian energy flow from the fc-th oscillator to the j-tli one, normalized with 
respect to A (see [BoLeR]). Or, for short, the energy flow. We prove 

6 Here and further on the rate of smallness of A is independent neither from e nor from the choice of the set C. 
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Theorem 1.6. Let p £,x be the unique stationary measure of eq. (1.6). Assume X to be 
sufficiently small. Then, under suitable assumptions for the interaction potential V, there 
exists a C l -smooth strictly positive function k : i —> M + satisfying n(x,y) = n(y,x), 

such that for every j,k G C, \j — k\ — 1, we have 

(//’ A , Jkj) —> Xnfll X, Tj)(Tk — Tj) + o(A) as £ —> 0 uniformly in C, (1.15) 

where o(A)/A —» 0 as A —> 0 uniformly in C. 

For details see Theorem 3.4. We emphasize that the function ft(7fc,7)), which we 
call the conductivity , depends only on the temperatures 71,7), and not on the choice of 
the set C. Convergence (1.15) shows that in the stationary regime the averaged flow of 
energy is approximately proportional (uniformly in C) to the local temperature gradient, 
if the noise and the interaction between oscillators are sufficiently weak. The coefficient 
of proportionality k is strictly positive and depends on the temperatures in a sufficiently 
smooth way. To prove Theorem 1.6 we note that Theorem 1.3 implies (/a £,x ,Jkj) —>■ 
(/i A , Jkj) as £ —» 0 uniformly in C, where /i A is the unique stationary measure of the 
effective equation (1.10). Then we analyse the development of the measure p x in A. 

Remark 1.7. For the reader interested in the limiting (as e —» 0) behaviour of the sta¬ 
tionary thermostatic energy flow, we note that it can be easily obtained from eq. (1.13). 
Indeed, let us average the both sides of (1.13) with respect to the stationary measure p £,x . 
The left-hand side vanishes, while the limiting behaviour of the first term in the right-hand 
side is given by Theorem 1.6. 

Now for simplicity we fix the dimension of the lattice d — 1 and put C = {0,1,..., IV}. 
Let us study the limit TV —> oo. For x G [0,1] denote j x := [xN] and J{x) := Jj x +ij x - 
If we put the oscillators to the interval [0,1], a j'-th oscillator to the point j/N, then the 
function J(x) will describe the energy flow between the nearest to the point x oscillators, 
so it can be seen as the energy flow through the point x. Choose a temperature profile 
Tj '■= T(j/N), where T(x) is a C 1 -smooth positive function, defined for x G [0,1]. Then 
Theorem 1.6 immediately implies 

lim hmhm^(//’ A , J{x)) = k{T{x))^-T{x), where k(y) := n{y,y). (1.16) 

iV—>ooA-^Oe—>0 A dX 

If the temperature profile is linear, i.e. T(x) = T\X + 7o(l — x), where 7o,7l >0, then 
(1.16) implies that 

lim lim lim -^-(//’ A , J(x)) — k(T(x)), (1.17) 

JV-s-oo A-U) £^0 XST 

where 5T := 7) — 7o- This resembles the Fourier law, see [BoLeR], 

Next we study stationary space-time correlations of the energy flow. In particular, we 
get the following result. 

Theorem 1.8. For any T > 0 and N gN we have 


1 T N -1 N -1 

k (f) = I™ J™ -fflj / (F £ ’\ J n+Wr X ( Jjj+i) > d T 

' iV n 1=0 1=0 


( 1 . 18 ) 


where *}3y A denotes the Markov semigroup associated with equation (1.6) with the temper¬ 
atures satisfying Tj = T for every j G C, while p £,x denotes its unique stationary measure 
(which is the Gibbs measure at the temperature T). 




See for details Section 3.3. Relation (1.18) resembles the Green-Kubo formula, see 
[Sp],[BoLeR] and [LepLiPo], ' To prove Theorem 1.8 we employ the averaging technique 
developed above, an explicit formula for the conductivity k, given in Theorem 3.4 (which 
is a full version of Theorem 1.6) and the uniformity in e of the rate of mixing for eq. (1.6). 

Despite that (1.17) and (1.18) resemble the Fourier law and the Green-Kubo formula, 
of course they are not. 

1.4 Assumptions and agreements 

Assumptions. Throughout the text we assume that the following hypotheses hold. 

HV. The function V : M 2 —> R is C r -smooth, r > 2, and V(x,y) = V(y,x) for any 
x,y G M. There exists a constant C > 0 such that 

\V(x,y)\<C(l + x 2 + y 2 ), (1.19) 

\d Zl V(x,y)\ <(7(1 + 1x1 + 12/1), \d 2 ZlZ2 V(x,y)\ <C(l + x 2 + y 2 ), (1.20) 

where Z\,Z 2 € {x,y}. Moreover, higher orders partial derivatives of V have at most a 
polynomial growth at infinity. 

HI. There exist constants «o > 0 and C > 0, independent from the choice of the set 
C, such that the initial conditions (p 0 , g 0 ) satisfy 

E e <*o(p 0 j+9<y) < (j f or every j £ C. 

Remark 1.9. We use relations (1.19)-(1.20) only to obtain uniform in C, e and r a-priori 
estimates (2.1) for solutions and stationary measures of eq. (1-6), and to prove that eq. 

(1.6) is mixing with independent from e rate. Once these properties are given, we do not 
more need (1.19)-(1.20). 

More specifically, for Theorem 1.1 we need the first estimate from (2.1) to be held. 
For Theorems 1.3 and 1.6 we need (2.1) and uniqueness of the stationary measure for eq. 

(1.6) . For Theorems l.f and 1.8 we need the same, and additionally, we need eq. (1.6) 
to be mixing with independent from e rate. 

Note that the exponential estimates (2.1) are superfluous. It suffices to have polynomial 
estimates of sufficiently large degree. Moreover, if we have only estimates which are not 
uniform in C then all the results remain true except the uniformity of convergences in C. 

Agreements. 1) We refer to item (i) of Theorem 1.4 as Theorem 1.47, etc. 

2) By (7, (7i, C* 2 ,... we denote various positive constants and by (7(a), (7i(a),... - 
positive constants which depend on a parameter a. Sometimes we skip the dependence 
of the constants on parameters. However, unless otherwise stated, we always indicate if 
they depend on the choice of the set C, times t, s,r ,..., positions j, k,l,m, ... e C and 
small parameters e, A. Constants (7, (7(a),... can change from formula to formula. 

3) We use notations a A b := min(a, b ), a V b = max(a, b). 

4) For MeN and vectors a = (a*,), b = (6 fc ) e C M , by a • b we denote the Euclidean 
scalar product in C M ~ 1R 2M , and by |a| the corresponding norm, 8 

a ■ b y^(Re a*, Re bk + Ima& Im bf) — Re afbk and |a| 2 := a • a = |afc| 2 - 

k k k 

7 In formula (1.18) we could also send N to infinity, however it does not make sense: limit (1.18) is the same for all N. 

8 We denote the / 2 -norm in the space of complex vectors C M and the Zi-norm in the space of indices Z d by the same sign 
| • |. This will not cause misunderstanding. 
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5) By saying that a function / : C M i —> R is C' n -smootli we mean that it is FA-smooth 
as a function / : M 2M i—> M. 

6 ) Let X be some vector space and / : ( Xj)j £ c £ X ^ i —> M. By Supp / we denote a 
minimal set S C C such that the function / is independent from (xj)j^s- 

7) For a metric space X by £b{X) (£i oc (X)) we denote the space of bounded Lipschitz 
continuous (locally Lipschitz continuous) functions from X to R. 

8 ) Assertions of the type "something is sufficiently small/sufficiently close to one/... " 
always assume an estimate independent from the choice of the set C. 

9) We assume 0 < £ < 1 to be sufficiently small where it is needed. 
Acknowledgments. I am very grateful to my PhD supervisors S. Kuksin and A. 

Shirikyan for formulation of the problem, guidance and encouragement. I would like 
to thank N. Cuneo, J.-P. Eckmann, V. Jaksic, J.L. Lebowitz and C. Liverani for useful 
discussions, and N. Cuneo and J.-P. Eckmann for the excellent stay in Geneva. This work 
is supported by the Russian Science Foundation under grant 14-50-00005 and performed 
in Steklov Mathematical Institute of Russian Academy of Science. 

2 The limiting dynamics 

In this section we study the limiting as £ —> 0 behaviour of solutions of equation (1.6). 
Most of results of this section are independent from the size of A, where 0 < A < 1. 

2.1 Estimates on solutions 

It is well known that eq. (1.6) has a unique solution and this solution is defined for all 
r > 0 (see [Khal2]). In Appendix A we prove that eq. (1.6) is mixing. First we need to 
obtain estimates for its solution and stationary measure, which are unifrom in e, r and C. 

Lemma 2.1. There exist constants a > 0,£o > 0 and C > 0 7 independent from the choice 
of the set C, such that for all 0 < e < Eq, 0 < A < 1, j 6 C and r > 0 we have 

(f)E max e °(p?«+9?M) < C, (ii) (p, e a ^ +q ^) < C, ( 2 . 1 ) 

se[r,r+l] 

where (p,g)(r) is a solution of eq. (1.6), satisfying V(p,q)( 0) = V(p 0 ,q 0 ), and p is its 
unique stationary measure. 

Proof. Item (ii) follows from item (i), the mixing property of eq. (1.6) and the Fatou 
lemma. 9 Let us now explain the idea of the proof of item (i). Fix some 1/2 < 7 < 1, 

Find nPDOTP 

U’(p, 9 )■■=’£, ( £ t(p, 9 ) + \pk 1 k) , j e z d , ( 2 . 2 ) 

fceC 

where we recall that the local energy £ k is defined in (1.12). 10 Estimate (2.1).i with 
the constant C depending on the number of nodes \C\ can be obtained in a standard 
way, by applying the Ito formula to the function e aU3 \ 1= \, see (2.9). To derive it, one 

should crucially use that for 7 = 1 the function := 7 ^ k ~^£k(p, q) coincides with the 

kec 

9 In fact, it can be obtained even without employing of the mixing property, in a slightly more complicated way (see 
Section 2.5.2 from [KS] ). One should use that dynamics "forgets" initial conditions (see (2.22)), and the Fatou lemma. 

10 For the brevity of notations here and further on in this proof we skip the upper index is. 
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Hamiltonain H, so that £ J | 7=1 is conserved by the Hamiltonian flow. The obtained in such 
a way estimate is not uniform in C because the r.h.s. of (2.9) contains terms proportional 

to the sum fff ^~ k \ which equals to \C\ if 7 = 1. To overcome this difficulty, we put 
kec 

7 < 1, so that the sum above is bounded uniformly in C. Despite that in this case the 
function is not a first integral of the Hamiltonain flow, the nearest-neighbour character 
of interaction between oscillators provides that if 7 ~ 1 , is its "approximate integral", 
in the sense that e~ 1 {£\ H} ~ (1 — 7 ), see (2.13). This suffices to establish the desired 
estimate. 

Now we start the proof. First we will establish the following auxiliary proposition. 
Introduce the family of norms 


n ■= 


y j 7 ^' k ^x\, where j e and x = (xk)k&c £ 


(2.3) 


fcec 


Proposition 2.2. There exists £7 > 0 such that for any 1/2 < 7 < 1, j 6 Z d and 
0 < £ < £1 we have 

(Ml? + IMI?)/4 < U’(p, q ) < C(||p||= + || 9 ||?) + <7,(7), ( M ) 6 R 2 ' c '. 

P 2 \ ^ 

Proof. Let us prove the upper bound. Since (1.19) implies £k + C V qf + C, 
we have 

Ui <Y,p- k 'fi + C V 9 2 + c + p|±«I) <^ 7 b-H ( p 2 + Ci V qf+C). 


k£C 


l:\l-k\<l 


k£C 


l:\l-k\<l 


Since 7 > 1/2, we have ^ 7 ^^ ^ qf <2^ ^ 7 < C||g||j. Moreover, 

k&C l:\l-k\<l k&C Z:|/-fc|<l 

Yf, < Y2 7^ < ( 7 ( 7 ). Then we obtain the desired estimate. The lower bound 

fcGC mG 

follows from (1.19) and smallness of e in a similar way. □ 

Now let us formulate the following result, which we will establish after the end of the 
proof of the lemma. 

Lemma 2.3. Let 1/2 < 7 < 1 be sufficiently close to one and 0 < £0 < £ i be sufficiently 
small, where £1 is defined in Proposition 2.2. Fix any <5 0 > 0. Then for 0 < 6 < ho 
sufficiently small and any solution (p,q)(r) of eq. (1.6), satisfying for some j G 7L d the 
estimate 

f)g { 0 W(p,q)(0) < £7 

we have 

E max ^ 

se[r,r+l] 

for any 0 < e < e 0 , 0 < A < 1 and r > 0, where C\ = C\{^,£q,8,C). 

We fix the parameters 7 and £ 0 as in Lemma 2.3 and do not indicate dependence of 
constants from them. Using first the upper bound from Proposition 2.2 and then the 
Jensen inequality joined with ass. HI, we see that 

BjW'K*] < c 1 E f ‘i c (l»l'tl«l?) < c 2 , (2.6) 


(2.4) 

(2.5) 


11 



for 5 0 sufficiently small and any j G 7L d . Since pj + q/ < \\p\\ 2 + ||g|| 2 , the lower bound 
from Proposition 2.2 joined with Lemma 2.3 implies that 

E max e a(p ? (s)+ ^ 2(s))/4 < E max e suHp ’ q){s) < C, (2.7) 

s6[t,t+1] s6[t,t+1] 

for 5 > 0 sufficiently small. Since the constant C 2 from (2.6) is independent from j, the 
constant C from (2.7) also is, so that we get the desired estimate with a := S/A. □ 

Proof of Lemma 2.3. Without loss of generality we assume that j = 0 and skip the 
index j in the notations for f/ J and || • ||j. Moreover, for simplicity we put A = 1: analysing 
the proof below it is easy to see that all estimates we obtain are uniform in 0 < A < 1. 

Step 1. First we will find 7,e 0 and > 0 such that for any 0 < £ < e 0 and s > 0 we 
have 

E e SlU ^<C, (2.8) 

where we have denoted U(s ) := U(p,q)(s). Further on we present only formal compu¬ 
tation which could be justified by standard stopping-time arguments (see, e.g., [KaSh]). 
Applying the Ito formula to e 5lU ^ we obtain 

H}-Y, P,d K e“ U( ‘ ] + £ + f -M 

zeC igC 

= s.e^iy, + y 2 + y 3 ) + ( 2 . 9 ) 

where 


y^s-^H}, y 2 = -J2pi d Pl l A ^3 = ^7I(a p 2 ? [/ + h 1 (a K f/) 2 ), 

zee i&c 


and the martingal M s ^ s := 5, / ^2Tid Pl Ud/3i, where so < s. Let us estimate 


so 


zee 


the terms 3d, 3 ^ 2 ,A 3 separately. 

Term 33- Writing 3d in more details we obtain 


n = ^ 1 E h] + e f-iPkik, h}. 


1*1 


fcec 


fceC 


We have 


e~'{£ k ,H}= {|,i E C(fc. 9 ,)} + {i E 


\l—m\=l 


l:\l—k\=l 


{y, E C(ft,9,)} + i E 


mGC 

p'k + Pi 


l:\l—k\=l 


l:\l—k\ = l 


E 

Z:|Z—fc|=l 


Pl-P 2 i 


,V(q k , 


( 2 . 10 ) 


( 2 . 11 ) 
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Note that for any ( x k i)k,iec e R |CxC| we have Y Y x ki = \ Y i x ki + x ik)- Using 


( 2 . 11 ) and applying the latter identity with x k i = 7 


fcec l:\l-k\=l p—fc|=l 

— — ,V(i we find 


= E E 7 W { 

keC l:\l-k\=l 


\k\\ Pk-Pl 


,V(q k , 


k£C 


5 E 


t m .yJ Pt-rf 


p_fc|=l 


.rfe,®)}- 


( 2 . 12 ) 


Since for k, l G C satisfying \l — k\ = 1 we have |/| — \k\ = ± 1 , we obtain 


| 7 l fc l_ 7 l*l| =7 lfc||i_ 7 l- 


“M: 1 , < 2771 - 7 ), 


7 |fc| 7 H 1 -7), if \l\ ~ \k\ = -1 
where we have used that 1/2 < 7 < 1. In view of estimate (1.20), we find 
, V(q kl <?;)) < C(p 2 k +p 2 i+q 2 k + qf + 1). Then from (2.12) we get 


Pi ~ Pi 


£ l ^i^{£k, H } <(i-i)c ^ kl (pl+p‘i+Qk+ ( i‘i+ 1 ) < ( 1 - 7 )c , i(lbll 2 +lkir)+c , 2 ( 7 )- 

fceC |i-fc|=i 

(2.13) 

7^ 

Let us estimate the term 2 . - {p k q k ,H} from (2.10). We have 

fcec 

{p k q k , H} = p k d Pk H - q k d qk H = p 2 k - q 2 k - £ ^ q k d qk V(q k , q t ). 

l:\l- k \ = l 

Then, using that \q k d qk V ( q k ,qi)\ < C(q 2 k + qf + 1), we hnd 

|fc| 1 

^2^2~{Pkqk,H}<-(\\p\\ 2 -\\q\\ 2 ) + Ce^2 7 |fc| (qI + Qi + 1) 

fcec fcec i-.\i-k\=i 

< |b|| 2 /2 + (-1/2 + Cie)\\q\\ 2 + 6 * 2 ( 7 ). (2-14) 

Combining (2.13) and (2.14) we obtain 

y < ((1 - 7 )C+ 1/2) INI 2 + ((1 - 7 + <0C - 1 / 2 ) 11 # + C,(7). (2.15) 

Term y^. Since 

9 pi U = (pi + 7^qi), (2-16) 

we hnd 

2 = -INI 2 - ^^i^piqi < (-i + NINI 2 + £|kll 2 - (2.17) 

A leC 

Term 3/a- In view of (2.16) we have d^U = 7 ^. Then, using that the temperatures Ti 
are bounded uniformly in l 6 C (see (1.2)), we hnd 

33s = Ti(^ w + <*i7 M {Pi + 3 ®) 2 ) ^ C (t) + 25i(||p|| 2 + Ikll 2 )- (2-18) 

lec 
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(2.19) 


Substituting (2.15),(2.17) and (2.18) to (2.9), we get 

2 < ^‘‘"(Adlpll 2 + ||,|| 2 ) + C(7)) + 4 m, 0 ,„ 

where 

A( 7 , e, hi) := ((1 — 'y'jC + 1/2 — 1 + £ + 2hi) V ((1 — 7 T £)C — 1/2 + £ + 2hi). 

Fix 1/2 < 7 < 1 sufficiently close to one, £ 0 > 0 and 0 < hi < ho sufficiently small, in such 
a way that A(y, e, hi) < 0 for £ < £q- Then, using Proposition 2.2, for e < £$ we obtain 

-^e SlU < hie 5 lC/ (A 0 U + C) + ^M SOiS , (2.20) 

for some A 0 < 0 , where we do not more indicate dependence of constants on the param¬ 
eters 7 ,£o and hi. Then for an appropriate constant C we get 

< _ e 5,u + c+ (2.21) 

as as 

Fixing So = 0, taking the expectation and applying the Gronwall inequality to (2.21), we 
obtain 

E e SlU(s) < E e 5lU{0) e~ s + C, (2.22) 

so that we have ( 2 . 8 ). 

Step 2. Now we will prove (2.5) for 0 < h < hi/2. Integrating inequality (2.21) with 
hi replaced by h over the interval [r, s], where t < s < r + 1 , putting So = r and taking 
E max , we get 

se[r,r+l] 


E max e 5u{s) <Ee su{r) +C + E max M r , s < C x + E max M r , s , (2.23) 

s6[t,t+ 1] se[r,r+l] ’ sS[r,r+l] 

where we have used (2.8). The Doob-Kolmogorov inequality ensures that 

E max M ts < CE J[M] TT+1 < cJ E [M] TiT+ i, (2.24) 

se[r,r+l] V V 

where [M] T)S denotes the quadratic variation of the martingale M rs . Due to (2.16), 

T+l T+l 

[M ] T , T+1 = + /e“^27I 7 2|,| (ft + f«) 2 *<C J e” £ 7 2 «(p 2 + ,?) is 


lec 


leC 


T+l 


< 4G / e 25u Uds , 


in view of the lower bound from Proposition 2.2 and the estimate ||p|| 2 +||g|| 2 > ^ 7 2 ^(p 2 + 

lec 

qf). Since 2h < hi, estimate (2.8) provides E [M] T)T+1 < C, so that (2.23) joined with 
(2.24) implies the desired relation (2.5). □ 
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2.2 a-variables and the effective equation 

Instead of studying eq. (1.6) in the real (p, ^-coordinates, it is convenient to pass to the 
complex variables 

u = ( u j)jec e C |c| , Uj(p, q ) = pj + iqj, j e C. 

Then the actions and angles take the form 

\ u f 

I = (Ij)jec, J j( u ) = and <P = (^j)jec, <Pj(u) = arg Uj. 

Denote 

V jk [u) := V(q j {u),q k {u)) = ^Imu^Imu*), 
where j,k G C satisfy \k — j\ = 1. Then in the u -variables eq. (1.6) takes the form 

iij = e^iuj + Pj(u) + \J2Tj/3j , (2.25) 

where 11 

P = ( Pj)jec, Pj{u) := A ^ WjVjkiu) ~ Reuj, (2.26) 

k£C:\j — k\ = l 

and Vj := 2d Sj is the gradient with respect to the Euclidean scalar product in C ~ M 2 . 
The derivative du : is understood formally, du j := (d Pj +id qj )/2. 

Let u £ (t ) be a unique solution of eq. (2.25) satisfying V(u £ (0)) = T>(u 0 ), where u 0 := 
Po + iQo- Our first aim is to study its limiting behaviour as £ —> 0. In the main order 
with respect to e eq. (2.25) describes a system of fast synchronously rotating uncoupled 
oscillators (see also </>equation from (1.8)). Let us pass to a rotating coordinate system 
where the oscillators are immovable. Put 

a £ = (a £ ) jeC , a £ (r) := e" l£_lr ^(r), j e C. (2.27) 

Obviously, we have 

\u £ (t)\ = |a £ (r)| for all j e C, and I{u s (t )) = J(a £ (r)). (2.28) 

The process a e (r) satishes 

a) = e~ ie ~ 1 ' T P j (e i£ ~ lr o £ ) + y/fr j e~ ie ~ lT $ j , j E C, (2.29) 

D(a £ (0 ))=D(n 0 ), 

where e l£ lr a £ is a vector with components (e l£ lr a e ) = e l£ lr a|, k G C. Regarding (2.29) 
as an equation in the extended phase space x R + 9 (a, 6), we have 

a £ = e~ ie Pj(e id a £ ) + y/Wje^Pj, 6 = e~ l . (2.30) 

Eq. (2.30) is a usual fast-slow system with the unique fast variable 6. Applying the 
method of stochastic averaging due to R.Khas’minski, we will show that its limiting (as 
£ — y 0) dynamics is given by the a-equation, averaged with respect to 6, see the next 
section. The latter has the form 

Uj = 77, (n) + y/7j(3j, j e C, (2.31) 

11 Note that the function Pj is real valued. 
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where 


2t r 


(2.32) 


7 ZAu) := [ e~ ie PAe ie u)ae , dO := —, 

.7 27T 

o 

and (3 = (f3f)j & c is a standard complex |C (-dimensional Brownian motion, so that (3j = 
P) + iff, where (P*) k =i, 2 , j^c are standard real independent Brownian motions. Following 
[KM], we call (2.31) the effective equation. 

Let us emphasize that under the limit e —> 0 the noise "bifurcates". Indeed, the rank 
of the dispersion matrix of the effective equation written in the real (p, g)-coordinates 
equals to 2|C| (so the diffusion matrix is non-degenerate), while that for eq. (2.29) equals 
to |C|. This happens because the fast rotation decomposes the noise into two independent 
components, see Proposition 4.3. 

Let us calculate the drift IZfu). For this purpose we will need to define the resonant 
averaging. Take a continuous function / : — y R and write it in the action-angle 

coordinates (I,(p). Then the resonant averaging (f)n is defined as the averaging of / in 
angles in the direction 1 = (1,1..., 1) of their fast rotation (see ^-equation in (1.8)). 
More precisely, 

2tt 

(. f)n{I,<p):= I f(I,<p + 61)d0. (2.33) 

o 

In the u- variables it takes the form 


2tt 

(f)R(u) = j f(e ie u)dO. 
0 


(2.34) 


Further on we will often use properties of the resonant averaging given in Appendix C. 
Denote by H res the resonant averaging of the interaction potential, 

HT “'=\ A V ^‘. where \r« := (2.35) 

\m—k\ = l 


We will call the function Vfff the resonant potential and H res — the resonant Hamiltonian. 
In view of (2.32) and (2.26), we have 


7 Zj ( u ) 



R e(e l6 Uj)dO 


(2.36) 


= A?'V, 


, 

- E 

2 ^ 

\m—k\=l 


2tt 


V mk (e lU u)d6 - e 


-i6 


iO i — 

e Wo + e Uj 


(10 = A/T,//' 


[u - -. 

v ; 2 


Then the effective equation (2.31) takes the form 

., A/V ; // r ' " - -y - v j 6 C, (2.37) 

and describes dynamics of a Hamiltonian system given by the resonant Hamiltonian XH res , 
where each particle is coupled with its own stochastic Langevin-type thermostat. Note 
that eq. (2.37), written in the real (p, ^-coordinates, coincides with eq. (1.10) from the 
introduction. 
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2.3 Averaging theorems 

It is known that eq. (2.31) has a unique solution, this solution is defined globally (see 
[Khal2]), and that eq. (2.31) is mixing (see [Khal2, Ver87, Ver97]) 12 . In Section 4.2 
we will prove the following theorem, which jointly with (2.28) immediately implies Theo¬ 
rem 1.1. 

Theorem 2.4. On the space C'([0, T], C^) we have the weak convergence of measures 

V(a e {-)) — 1 V(u{-)) as e —» 0 uniformly in C, (2.38) 

where u(t) is a unique solution of the effective equation (2.31 ) satisfying T>(u(0)) = T>{uq). 

Recall that the uniformity in C of the weak convergences of measures through all 
the text is understood in the sense of finite-dimensional projections. For example, for 
convergence (2.38) it means that for any bounded set A C 7L d and any continuous bounded 
functional / : C([0, T], (d c l) i->- M, satisfying 13 Supp / C A, we have 

E/(a e (-)) —» E/(«(•)) as £ —> 0 uniformly in C, 

where we assume that the set C satisfies C D A. 

Convergence (2.38) transfers the estimates for solutions of eq. (1-6) obtained in 
Lemma 2.1 to the effective equation. 

Proposition 2.5. Let u(r) be the unique solution of the effective equation ( 2.31 ) satisfying 
V{u{ff)) = V[uq), and /i be its unique stationary measure. Then for all 0 < A < 1, j G C 
and t > 0 we have 

(i) E max e a \ u ii s )\ 2 < C , {ii) (//, e “ |u < C, 

sS[r,r+l] 

where the constants C and a are the same that in Lemma 2.1. 

Proof. Let u £ (t ) be the unique solution of eq. (2.25) satisfying D{u £ ( 0)) = V{uq). Writ¬ 
ing the first estimate from Lemma 2.1 in the complex variables we get E max e“^F)l < 

s€[r,r+l] 

C, for all 0 < £ < Eq, 0 < A < 1, j £ C and r > 0. Then (2.28) joined with Theorem 2.4 
and the Fatou lemma implies item (i). Item (ii) follows from item (i), the mixing property 
of the effective equation and the Fatou lemma. □ 

Next we investigate the limiting (as £ —y 0) behaviour of the unique stationary measure 
p £ of eq. (2.25). Let us introduce an additional assumption. Consider an unbounded set 
Coo Q Tf. Take the inductive limit as C Coo of the effective equation, i.e. eq. (2.31) 
with the set C replaced by Coo■ We will call it the Coo-effective equation. 

Definition 2.6. (i) A weak solution u°°(r) of the Coo-effective equation is called regular 
if there exists A > 0 such that for any T > 0 we have 

sup E max e A l u rC)l 2 < qq. 

jeCoo o < t<t 

(ii) A probability measure p°° defined on the Borel u-algebra B(0 Co °\) 14 is called regular 
stationary measure of the Coo-effective equation if the latter admits a regular weak solution 
u°°(t) satisfying V{u°°{t)) = /i°°. 

12 In order to establish the mixing property one could also use the method explained in Appendix A. 

13 Recall that Supp / is defined in Agreements. 6. 

14 The space Cl c °°l is provided with the Tikhonov topology of pointwise convergence. 
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Sometimes we will assume 

HCqo. For each unbounded set C 7L d the corresponding C^-effective equation admits 
a unique regular stationary measure. 

Passing to the limit C C^ and using estimates of Proposition 2.5, we can prove that 
the Coo-effective equation has a regular stationary measure, but we do not know if it is 
unique (however, we believe that it is). The following result provides a sufficient condition 
for fulfillment of ass. HCoo. 

Proposition 2.7. Assume that A is sufficiently small and second partial derivatives of 
the resonant potentials Vff s are bounded. 15 Then assumption HCoo is satisfied. 

In particular, if the second partial derivatives of the interaction potential V are bounded 
then those of Vff s also are. Proposition 2.7 is a corollary of Proposition 4.9, where we 
show that under conditions of Proposition 2.7 the Coo-effective equation defines a mixing 
Markov process, in some appropriate space. Its proof is based on a simple observation 
that for A sufficiently small the Coo-effective equation is contracting there. See Section 4.4. 

Theorem 2.8. (i) Let ff and /i be the unique stationary measures of equation (2.25) and 
the effective equation (2.31) correspondingly. Then // —^ p, as £ —> 0. 

(ii) If assumption HCoo Is satisfied then the convergence above holds uniformly in C. 

The proof of the theorem is postponed to Section 4.3. It is based on a version of 
Theorem 2.4 for stationary solutions, and rotation invariance of the effective equation 
discussed below. 

In Appendix A it is proven that the rate of mixing of eq. (2.25) is independent from 
e. Jointly with Theorems 1.1 and 2.8, this implies the following result. 

Theorem 2.9. (i) The convergence T>(l(u £ (r ))) —^ T>(l(u(r))) as e —» 0 provided by 
Theorem 1.1 holds uniformly in r > 0 in the sense that for any continuous bounded 
function f : i —y M we have 

sup |E/(J e (r))-E/(/(r))| ^0 as £ -> 0, (2.39) 

T> 0 

where I £ (t) := I(u £ (r )) and J(r) := I(u(r)). 16 
(ii) We have 

lim lim V(I £ (t)) = lim limT , (/ e (r)) = lim T>(l £ (r)) = II/*/i, 

£—>•0 T— >• OO T— >• OO £—^0 £—>-0 

T —^OC 

where p is the unique stationary measure of the effective equation (2.31), and II/ denotes 
the projection to the space of actions. 

Proof, (i) Eq. (2.25) and (2.31) are mixing and the rate of mixing for eq. (2.25) is 
uniform in e < £i (see Theorem A.l). Then for any <5 > 0 there exists s > 0 such that for 
any £ satisfying £ < eq 

sup |E/(/ £ (t)) — (II/*//, /)| < S/3, sup|E/(/(r)) -(n/^,/>| <5/3, (2.40) 

T>S T>S 

15 If the second partial derivatives of the function V£j 8 are bounded for one pair of indices k,j E C, then they are bounded 
by the same constant for all other k,j E C. Indeed, this is evident since V£j S (u) = (V(qk( u )i Qj( u )))R- 
16 We do not know if the convergence holds uniformly in both C and r. 
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where fj£ is the unique stationary measure of eq. (2.25). In view of Theorem 2.8 .i, there 
exists £2 such that for any e < e 2 we have 

|(n 7 *//,/) - (n 7 */q/)| < 5/3. (2.41) 

Inequalities (2.40) and (2.41) imply that for £ < £1 A £2 we have 

sup |E/(/ e (r)) — E/(/(t)) I <5. (2.42) 

T>S 

Applying Theorem 1.1 and using that the weak convergence of measures is equivalent to 
convergence in the dual-Lischitz norm (see [Dud], Theorem 11.3.3), it is not difficult to 
show that sup | E / (l £ (t))-E/(/(t))|-> 0 as £ —> 0. Then there exists £ 3 , such that 

0<T<S 

for any £ < £3 we have sup IE — E / (/( t)) I < 8. Jointly with (2.42) this implies 

0 <T<S 

that for £ < £1 A £2 A £3 we have sup |E /(/ e (r)) — E/(/(r)) I < 5. 

T> 0 

(ii) The hrst two convergences follow from the mixing properties of eq. (2.25) and 
(2.31), and Theorems 1.1,2. 8 .*. The last convergence follows from item fz) joined with the 
mixing property of eq. (2.31). □ 

2.4 The limiting behaviour of the local energy 

Recall that the local energy £f(p, q) is defined in (1.12). We will consider it as a function 
of the complex variables, and abusing notations we will write £'■ (u). Define the vector of 
local energy 

£"(«) := 

Let u £ (t) be a solution of eq. (2.25), either satisfying T>(u £ ( 0)) = T>(uf) or a stationary 
one. Since £ u = / + O(u), the limiting as £ —» 0 behaviour of the vector of local energy 
£ v (u £ ) coincides with that of the vector of actions I(u £ ), while the latter is governed by 
the effective equation, in the sense of Theorems 1.1,2.8 and 2.9. More precisely, we prove 

Proposition 2.10. Letu £ (r) be a solution of eq. (2.25) satisfying T>(u(0)) = T>(uf), and 
u £ (t) be its stationary solution. Let u(t) be a solution of the effective equation (2.31) 
satisfying T>(uf ))) = T>(uf), and /a be its unique stationary measure. Then, denoting 
£ u (t) := £ v (u £ {t )), £ u (t) := £ u (u £ (t)) and J(r) := I(u(t)), we have 

(i) T>{£ u f)) —*■!)(/(•)) as e—>■ 0 on C([0, T], RID) uniformly in C. 

(ii) V(£ u {t)) £>(/( t)) as £ —» 0 uniformly in r > 0 in the sense explained in Theo¬ 

rem 2.9.~ 

(Hi) T>(£ v (t)) —*■ Hup. as £ — > 0 for any r > 0. If ass. H0 o 0 is satisfied, this convergence 
holds uniformly in C. 

(iv) lim lim V(£ v (t)) = lim lim T>(£ v (r)) = lim V(£ v (r)) = IIj*/i. 

£—^0 t—»■ 00 v 7 r—>■ 00 e^-0 v 7 s —^0 v 7 

T—> OO 

Proof. Since the interaction potential V has at most a polynomial growth, the estimate 
(i) of Lemma 2.1 implies that the limit as £ —> 0 of T>(£ v ()) coincides with that of 
V(I{ u £ ())) in the sense that if one limit exists then another one exists as well and the 
two are equal. Moreover, if one convergence is uniform in C then another one also is. 
That is why item (%) follows from Theorem 1.1. Items (ii)-(iv) can be obtained similarly, 
using Theorems 2.9.Z, 2.8 and 2.9. ii correspondingly. For the hrst limit from item (iv) one 
should additionally use that eq. (2.25) is mixing. □ 
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2.5 Rotation invariance of the effective equation 

Further on we will need the following results. Since the effective equation (2.31) is obtained 
by the resonant averaging procedure, it is rotation invariant in the direction 1 of the fast 
rotation of angles: 

Proposition 2.11. Ifu(r) is a weak solution of the effective equation (2.31) then for any 
£ G [0, 2n) the process e^u(r) is also its weak solution. 

Proof. Process w{r ) = e^u(r) is a weak solution of equation 

1-K 

Wj = f e - l(0 ^ ) P J (e ,/1 u)d9 + v 7jr'%. j e C. 

o 

Changing the variable of integration 6 := 6 — £ and noting that {e l ^(3j, j G C } is another 
set of standard complex independent Brownian motions, we obtain that the process w(t) 
is a weak solution of the effective equation (2.31). □ 

Proposition 2.12. Let fi be the unique stationary measure of the effective equation (2.31 ). 
Then for any continuous function f : i —> M with at most a polynomial growth at infinity 

we have 

(dJ) = (d, (Dis¬ 
proof. Since the stationary measure /i of the effective equation is unique, Proposi¬ 
tion 2.11 implies that it is rotation invariant in the direction 1: for any £ G [0,27t) we 
have (e*^)*/i = fi. Then 

27T 27T 

(p, f) = j /}$; = (/b J f(e^-)df) = </i, (f) R ). 

o o 

□ 

3 The energy transport 

In this section we assume that A 1 and investigate the limiting (as £ —» 0) behaviour of 
the stationary Hamiltonian energy flow for eq. (2.25). We prove Theorems 3.4 and 3.12 
which are full versions of Theorems 1.6 and 1.8 given in the introduction. 

3.1 The limiting behaviour of the energy flow 

Consider the Ornstein-Uhlenbeck process (2.3T) | a=o 

u m = ——f 'J'Tm.ftm-, m G C. (3-1) 

It is well known that the unique stationary measure pP of eq. (3.1) has the form 

// ° = e“ |u| rd m , (3.2) 
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where \u\?j- := (27m) 1 \u m \ 2 and dm is a normalized Lebesgue measure on ~ M 2 I C I, 

mgC 

dm = dLeb/((27r)l t In me c7^„) • The Ornstein-Uhlenbeck process (3.1) is important for our 
study because of the following 

Proposition 3.1. (i) Letu x {r) and u(t) be solutions of the effective equation (2.37) and 
eq. (3.1) correspondingly, satisfying V{u x {fS)) = T>(u( 0)) =V(uq). Then 

V(u x (-)) —>■ V(u(-)) as A —^ 0 on C([0, T], C^) uniformly in C. 

(ii) Let p x be the unique stationary measure of eq. (2.37). Then p x —^ /i° as A —> 0 
uniformly in C. 

Proof. The required convergences can be established by arguments similar to those 
used in the proofs of Theorems 2.4 and 2.8, but significantly simplified. The uniformity 
of convergence from item (ii), similarly to that of Theorem 2.8. ii, follows from fulfillment 
of assumption analogous to HCoo. Namely, since equations (3.1) are diagonal, for any 
unbounded set Coo C Z d a stationary measure of their inductive limit as C /* (given 
by (3.1) with C = Coo ) is unique. □ 

For k, j G if denote by A k i the generator of the process (3.1) with C replaced by the 
set {k, j}. Considering A kj as an operator on real valued functions, we have 

A ti =Y,( 27?^ - | • V,). (3.3) 

1=3,k 

We will need the following well known property of the Ornstein-Uhlenbeck process. 

Lemma 3.2. (i) Fix k,j G C. Let a measurable function £ : i —> R. has at most a 

polynomial growth at infinity, satisfies Supp£ C {k,j} and = 0. Then equation 

A kj V = £ (3.4) 

has a unique solution rj : i-> M in a class of functions satisfying the following properties 

(a), ( b ), (c): (a) g belongs to each Sobolev class W 2 loc , p > 1, (b) r) has at most a polynomial 
growth at infinity, and (c) (y u°,r /) = 0. The solution g satisfies Suppry C {j,k} and has 
the form 

OO 

g(u) = - J ¥° t £(u) dr, (3.5) 

o 

where ^ denotes a Markov semigroup associated with the Ornstein-Uhlenbeck process 
(3.1). Moreover, the gradient \7rj(u) also has at most a polynomial growth at infinity. 

(ii) If the function £ is C n -smooth for some n G N, then the solution g is also C n - 
smooth. 

Proof, (i) This is a particular case of Theorem 1 from [PV] . 

(n) It follows in a standard way from smoothness of coefficients of the generator A k \ 
C' n -smoothness of the function £ and the Sobolev embedding theorem (see [Tay], Theo¬ 
rem 5.11.1). □ 

Recall that the energy flow Jkj{p,q) is defined in (1.14). Abusing notations, we put 
Jkj(u) := Jkj(p(u),q(u)). Introduce the resonant energy flow as 

Jtj‘ ■= (A;)*. 
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Proposition C.l .ii joined with ass. HV implies that the function Jf*f s is C r 1 -smooth, 
and all its partial derivatives have at most a polynomial growth at infinity. 

Proposition 3.3. We have 

(0 JIT = 2 d n v^ = - 2 d vj vg s - (ii) </x°, J™) = 0 . 

Proof, (i) The Poisson bracket written in the action angle coordinates (/, p) takes the 

form {hi,h 2 } = Yj{^ k hidi k h 2 - d Ik hid Vk h 2 ). Then (1.14) implies 
kec 


Jkj — {Ij Ik,Vkj} — (d Vk Vkj d Vj Vkj). (3-6) 

Proposition C.l .Hi provides that for l G {j, k) we have ( d Vl Vkj)R = <9^(14 j)n = 

It remains to take the resonant averaging of (3.6) and note that Proposition C.l .iv implies 


= -wr- w 

(ii) It follows from item (i) and rotation invariance of the Gaussian measure /i°. □ 

In view of item (ii) of Proposition 3.3, the resonant energy flow satisfies conditions 
of Lemma 3.2. Denote by rjkj(u ) obtained there solution of eq. 

A ki -n ki = JZT- ( 3 . 8 ) 


Theorem 3.4. Let p e,x be the unique stationary measure of eq. (2.25). Then for every 
j, k G C, | j — k | = 1, we have 

(/4’ A , ffkj) —^ Aft(7fe,7j)(7fc — Tj) + o(A) as e —> 0, (3-9) 


where o( A)/A —> 0 as A —> 0 uniformly in C, and 


K(T k ,Tj) 


_ (fi^TklVkVktf + Tj\Vm\ 2 ) 

T k Tj 2T k Tj 


(3.10) 


If assumption HCis satisfied then convergence (3.9) holds uniformly in C. 

The function k : R+ K > R satisfies n(x,y) = n(y,x), is C r ~ 1 -smooth and nonnegative. 
If the resonant potential Vjff s is not independent from the angles p, then k is strictly 
positive. 


Emphasize that the function n(Tk,Tj), which we call the conductivity , implicitly de¬ 
pends on the temperatures through the measure /i° and the function r/kj (the latter de¬ 
pends on them since the operator Al fcj does). 

Since rj k j is the solution of eq. (3.8) obtained in Lemma 3.2, formula (3.5) joined with 
(3.10) implies that the conductivity k can be also represented in the form 


<%,%) = J ^jg’dr). 

3 o 


Note that the sufficient condition for fulfillment of ass. HC 00 provided by Proposition 2.7 
is relevant for Theorem 3.4 since convergence (3.9) gives some information only if A is 

small. 

In order to prove Theorem 3.4 we will need the following proposition, which can be 
obtained by simple integration by parts. 
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Proposition 3.5. Let functions rf-^rf e iP|; oc (Cl c l,M) with their gradients Vt/,V? 7 2 
have at most a polynomial growth at infinity. Then 

(n° t r) 2 A ki r)') = -i (/, V TiV„i' ■ V,r f 2 > = 

l=k,j 


Proof of Theorem 3.4- Denote by the unique stationary measure of the effective 
equation (2.37). Theorem 2.8. * implies 

(/T' A , Jkj) ~> = (V‘ X >J'kj S ) as £ —^ 0 , (3-11) 

where the last equality follows from Proposition 2.12. If ass. HCoo is satisfied then 
Theorem 2.8 . ii implies that convergence (3.11) holds uniformly in C. Let us show that 

(T\Jk? s ) = A«(7fc,7))(7fe — 7j) + o(A), (3.12) 

where o(A)/A —> 0 as A —>■ 0 uniformly in C. The key role will play the following 

Lemma 3.6. Let the function £ e C 1 (Cl c l,M) satisfies Supp £ C {k : j}, has at most a 
polynomial growth at infinity and (//*,£) = 0. Then 

</i\£) = A</i 0 ,/£> + o(A), (3.13) 


where o( A)/A —> 0 as A —> 0 uniformly in C, f = \ fmn, where f mn = / nm; and 


fmn is the solution of eg. 


obtained in Lemma 3.2. 17 


\m—n\ = l 


T -T 

Amur _ 'li • m rrres 

S* J mn — ry ry J mn 1 

/ m / n 


(3.14) 


Before proving the lemma we will finish the proof of the theorem. Applying Lemma 3.6 
to the function £ = J[f s we get 

= \ E (3.15) 

\m—n\=l 

Proposition 3.7. If the sets {m, n} and {k,j} do not coincide then 

(AfmnJ£?) = 0. (3.16) 

Proof. Since Supp f mn C {m, n}, if m, n ^ k or m, n ^ j then the function f mn does not 
depend on the angle ip k or ipj correspondingly. Then (3.16) follows from Proposition 3.3 . i 
and rotation invariance of the Gaussian measure /i°. □ 

Since fjk = fkj, (3.15) joined with Proposition 3.7 implies 

</A J™) = A(/i°, f kj J™) + o( A). (3.17) 

Let rjkj be the solution of equation (3.8), obtained in Lemma 3.2. Eq. (3.14) implies 

fkj = W so that (/x°, fkjJfj 3 ) = k kj (T)(T k - 7)), (3.18) 

/fc lj 

17 Fix a bounded set A C Zf Almost without changes in the proof, a similar result can be established for a function £ 
satisfying Supp £ C A. We restrict ourselves to the case A = {k,j} since we will apply the lemma only for a two-dimensional 
situation. 
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where by T we have denoted the temperature profile (' Tj)jeC , and k k j(T) ■ = 


7771 


Using that the function Jff s does not depend on the temperatures T, rj k j depends on them 
only through the components Tk,Tj, that Suppf7^/ S ; Supp rj k j C {j,k} and pfi is a prod¬ 
uct measure, we see that the function k k j depends on T only through Tk,Tj, so that 
kkj{T ) = Kkj(7k, 7 j), for some function K k j. Obviously, the latter does not depend on the 
choice of the indices k,j, so we skip them and get 


k ( 77 , Tj) = 


Jtp 

7771 


(3.19) 


Now (3.18) joined with (3.17) and (3.11) implies the desired convergence (3.9). To obtain 
the second representation from (3.10) for the conductivity k, first we use eq. (3.8) and 
then Proposition 3.5: 


(n°,ri k jJ%‘) = W.miAVrto) = -iV Ti\v m \ 2 ), 

l=k,j 


(to, E T,\v lVtj \ 2 ) 

l=k,j 

so that we get n(T k ,Tj) = - —— -. In particular, this implies that the con- 

277 T) 

ductivity k is nonnegative. Moreover, if the resonant potential Vff s is not independent 
from the angles 99 , then in view of Proposition 3.3. * the function does not vanish. 
Consequently, r] k j is not a constant, so k is strictly positive for any Tk, Tj > 0. 

Let us show that n(Tk,Tj) = n{Tj,T k ). Since Jf.f s = —JJ k s and A k i = A ]k . we have 
Vkj = -Vjk • Then, due to (3.19), K(T k ,Tj) = -(777}) - 1 (iu 0 , VjkJJk 8 ) = «( 7j,T k )• 

It remains to prove that the function k is C' r ^ 1 -smooth. It is convenient to change the 
variables: let Vj = Tj l ^ 2 Uj, j £ C. Then, due to (3.19), we have 

K (7k,Tj) = -iTkTjY 1 (Jl 0 , J k j S rjkj), 

where /7° := e -M 2 / 2|c| dm, dm := dLeb/(27r)l c l, Jff s (v) := Jff s (u(v )) and rj k j{v ) : = 
VkjHv)). 

Proposition 3.8. For every v £ Cl c l the functions Jff s {y) = Jff s {v, T k , Tj) and rj k j(v ) = 
r)k](v, %, Tf) are C r ~ 1 -smooth with respect to the temperatures (77, Tj) £ M+. The corre¬ 
sponding partial derivatives till the order r — 1 inclusive are continuous in (v, 77, Tj) £ 
C' c l x and have at most a polynomial growth at infinity in v, which is uniform with 
respect to the temperatures ( 77 , Tj) if they belong to a compact set in K+. 

Since the function is independent from the temperatures 77, Tj, assertion of Propo¬ 
sition 3.8 concerning the function Jf ( f s is obvious. That, concerning rj k j is not so clear. 
It is not implied by standard arguments of elliptic theory because the generator A k i has 
growing coefficients. Its proof is based on simple analysis of formula (3.5) and is postponed 
to the end of the section. 

Since the measure JjP has an exponentially decaying density, Proposition 3.8 implies 
that the function (ft 0 , J k f s fikj) is C 7 " - 1 -smooth, so that the conductivity k also is. □ 
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Proof of Lemma 3.6. It is known (see [BKR]) that the unique stationary measure /i A 
of the effective equation (2.37) has a continuous density p x with respect to the normalized 
Lebesgue measure dm, so that dp x = p x dm. Then, in view of (3.2), we have 

dp x = h x dp°, (3.20) 

where the density h x satisfies 

p x = e -\n\ f h x ( 3 . 21 ) 

The proof of the lemma is based on a decomposition of the density h x in the small 
parameter A. Let us start with writing down an equation for h x . Denote by A the 
generator of the Ornstein-Uhlenbeck process (3.1), 

A = A T -fv, where A T := 2 £ T m ^ mE „. 18 (3.22) 

mec 

Then the generator of the effective equation (2.37) takes the form A + XB, where 
B:=g-V and g = (g m ) mGC , g m := iV m H res = i ^ 

n:\n-m\=l 

The density p x satishes in the weak sense the inverse Kolmogorov equation 

(A* + XB*)p x = 0, (3.23) 

where 

A* := A r + ^ ■ V + |C| and B* = -g-V. 

In order to calculate the operator B*, we have used the equality V • g — 0, following from 
the identity V m • = 0, which holds for any function tf G C 2 (Cl c l,M). Applying 

the operator A* to the both sides of (3.21), by the direct computation we find A*p x = 
g-Mr Ah x . Then, substituting (3.21) in (3.23) and multiplying the resulting equation by 
e Mr we get 

(A + A B*)h x = 0, where B*h x := e l <B*(e~ lul rh x ). (3.24) 

Now we will decompose equation (3.24) in the small parameter A, then find a solution 
of the obtained system till order A and show that the rest is small. Let us start with a 
formal computation. Put 

h x = 1 + A/ + A J A , (3.25) 

where / does not depend on A and f x = o(l) as A —* 0. Substituting (3.25) in (3.24), in 
order A 0 we have the identity A 1 = 0, while in the order A we obtain 

Af = -5*1. (3.26) 

To find the function fx, we solve eq. (3.26) and put 

/ A := A _1 (/i A - 1 - A/), (3.27) 

so that (3.25) becomes an identity. 

18 In the real ( p, g)-coordinates the operator A 7 ” has the form of the weighted laplacian A 7 " = 5 57 Tm(d 2 2 + d 2 2 )• 

meC Pm 
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Proposition 3.9. (%) Eq. (3.26) has a solution f of the form claimed in the formulation 
of the lemma. 

(ii) The function f x , defined in (3.27) with f from item (i), satisfies 

(A, / A £) —•► 0 as A —y 0 uniformly in C, (3.28) 


where £ is a function satisfying the conditions of the lemma. 

Take / and f x as in Proposition 3.9. Then, using that (A,£) = 0, (3.20) and (3.25), 
we obtain 

(A, £> = (A, h x 0 = A (/A /£) + o( A), (3.29) 

where o(A) = X(/jP, / A £), so that we get the desired decomposition (3.13). □ 

Proof of Proposition 3.9. Item (i). Using that for any function if(u) G C' 1 (Cl c l,M) we 
have (i'Vmfj) ■ u m = —AmA and the symmetry (3.7), we calculate 

= = C’%=-2 E 

m£C m£C n:\n-m\=l \n—m\ = l 

= E CC 1 - = \ E f rp-T- 1 )^, (3.30) 

\n—m\=l \n—m\=l 


where in the last equality we have employed Proposition 3.3. i. Take the solution f mn of 


a mn £ _ ( / ~T~ 1 _ ' 7 “—1 \ sjres _ 

71 J 17171 — \ ' m ’ n J'Jmn ~ 


Tn-Tm 

T n Tm 


V 

rr 


(3.31) 


obtained in Lemma 3.2. Put 

f ■ = 2 ^ ^ fmn• 

\m—n\ = l 

Since Supp f mn C {m,n}, we obviously have Af mn = A mn f mn . Then (3.31) joined with 
(3.30) implies (3.26). Since Jfff = —Jfff and A mn = A nm , equations for f mn and f nrn 
coincide, so that we get f mn = f nm . 

Item (ii). We will need the following two identities. Take any function 77 G fU 2 2 Zoc (Cl c l, M) 
such that 77 and V/7 have at most a polynomial growth. 

1) Since /i A is a stationary measure of the effective equation (2.37) and A + A B is the 
generator of the latter, we have ( A + A 13)*gf = 0. So that, (/A At 7) = — A(/i A , Bp). 

2) The integration by parts implies (p°,B*lp) = (/a 0 , Bp). 

Also let us note that Proposition 3.5 still holds if replace the operator A kj by A and 
the sum over l G {k,j} by the sum over l G C. 

Now let us take the solution p of equation A k ^p = £, obtained in Lemma 3.2. Since 
Supp 77 C {k,j}, we have A?7 = A k] p = f. Using (3.27), (3.20), the version of Proposition 
3.5 above, (3.26) and the identities 1), 2) above, we get 


(/A /A) = (A, f X Ap) = A h x Ap) - (A, (1 + Xf)Ap)) 

= a -1 ((AAA - <AA(1 + a/)A) = ~(h x ,Bg) + (fi°,B*ip) = -(A - AAA- 

Since the function £ is A-smooth, Lemma 3.2. ii implies that the solution p also is. So 
that, the function Bp is continuous. Then Proposition 3.1 .ii implies 

(A - A, Bp) 0 as A ->• 0. (3.32) 
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Recall that the uniformity in C of convergence in Proposition 3.1 .ii is understood, as usual, 
in the sense of finite-dimensional projections. Then, in order to prove that convergence 
(3.32) holds uniformly in C, it suffices to find a bounded set A C independent from 
the choice of the set C, such that Supp Brj C A. Since Supp ?7 C {k,j}, we have Brj = 
9j ' ^jV + 9k ' Vfc? 7 . Since for any m we have Supp g rn C {n : \m — n\ < 1}, the set A can 
be chosen as A = {n G Z d : \n — j\ A \n — k\ < 1}. □ 

Proof of Proposition 3.8. In this proof we accept the following agreement: saying that 
something is "uniform in T", we mean the uniformity in the temperatures ( Tk , 7j), if they 
belong to a compact subset of . 

Since the resonant energy flow is independent from the temperatures, Jf ( f s de¬ 
pends on them only through the function u(v), so the desired properties of the function 
Jff s follow. Let us turn to the function rjkj- Since r]kj is a solution of eq. (3.8) obtained 
in Lemma 3.2, due to formula (3.5) we have 

OO 

>>*(“) = - ( 3 . 33 ) 

0 

Since the change of variables transforms the Ornstein-Uhlenbeck process (3.1) to 

Vj = — A + $j, j G C, (3.34) 

relation (3.33) implies 

OO 

ij tl (v)=-J^ T jg‘(v)dT, (3.35) 

0 

where is a Markov semigroup associated with the process (3.34). Note that the 
operator ^ does not depend on the temperatures. Put dj- := d ll ^ l2 l2 , where 0 < l\ + I 2 < 

k j 

r — 1 . To prove that the function rjkj is (r — l)-time differentiable in Tk, Tj, we will just 
show that 

OO 

drVkj(v ) = - J ^d r Jff s (v) dr , (3.36) 

0 

where the last integral converges. Changing the variables u 1 —>■ v in the identity of Propo¬ 
sition 3.3 .ii, we get = 0. This holds for any Tk,Tj > 0, so that we have 

(p°, drJff s ) = 0. Since /i° is a unique stationary measure of the Ornstein-Uhlenbeck 
process (3.1), the measure Ji° is that for eq. (3.34). The mixing property of the latter 
implies 

|USrJfr(0| = |- (TP,d T J T t ‘‘)\< C( 1 + M”)e-y (3.37) 

for any r > 0 and some constants C,m,b > 0. It is possible to show that the constants 
C, m and b may be chosen uniformly in T. It follows that the integral from the r.h.s. of 

OO__ 

(3.36) converges uniformly in T. Then it is not difficult, to show that f Jlf s dr = 

0 

OO ^ _ 

dr J So that, applying dr to the both sides of (3.35), we get (3.36). 

0 
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Relation (3.36) joined with (3.37) implies that the function dq-rj k j{v) has at most a 
polynomial growth in v, which is uniform in T. The proof of the fact that the function 
dq-rjkj is continuous with respect to (v,Tic,Tj) is not complicated. Using that solutions of 
eq. (3.34) depend on initial conditions in a continuous way, one should show that the r.li.s. 

oo T 

of (3.36) with the time integral f ... dr replaced by f ... dr is continuous in (w, T k , Tj), 

o o 

OO 

and then observe that the rest f ... dr is small, in view of (3.37). □ 

T 


3.2 Examples 

1. Consider the quadratic interaction potential 

V(qj,q k ) = (qj-qk) 2 - 

It turns out that in this case the conductivity k(7), Tk) from Theorem 3.4 does not depend 
on the temperatures, so is a positive constant. We skip the proof since next we will consider 
a similar but more complicated situation. 

2. Let us now calculate the conductivity for the interaction potential 

v(qj,qk) = (qj -q k ) 4 - (3.38) 

The potential V has growth of power four, so assumption HV is violated. Despite this we 
have the following result. 

Proposition 3.10. (i) Eq. (2.25) is exponentially mixing with uniform in e rate (in the 
sense as in Theorem A.l). 

(ii) There exists e 0 such that for any 0<£<e 0 ;0<A<l ; r>0 and m > 0 we have 
E max (|p(s)| 2 + |g(s)| 2 ) m < C(\C\,m), (/i, (\p\ 2 + Vl\ 2 ) m ) < C{\C\, m), 

se[r,r+l] 

where (p,q)(r) is a solution of eq. (2.25) satisfying V(p,q)( 0) = T>(p 0 ,q 0 ), and p, is its 
unique stationary measure. 

Proof, (i) By the same arguments that was used in the proof of Theorem A.l. 

(ii) By applying the Ito formula to the function ( H u + | Pjqf) m ■ 19 □ 

j&c 

As it was explained in Remark 1.9, Proposition 3.10 implies that all theorems proved 
in this and previous chapters remain valid, except the uniformity of convergences in C 
which fails. 

Proposition 3.11. The conductivity k has the form 

K(T J ,T k ) = C(T i + T k ) 2 , 

where the constant C is independent from the temperatures. 

Proof. By the direct computation we get 

VS“(«) = ||«j—«*| 4 . 

19 For the both items one should use that (i ^2 P:i ( lj ■ § T, ( ( lk ~ 9j) 4 } = ~ E (ik ~~ 9j) 4 - 

jec |fc-i|=l |fc-il=l 
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To calculate the conductivity k, defined by formula (3.10), we need to find the solution 
r]kj of eq. (3.8), satisfying properties (a)-(c) from Lemma 3.2. Writing the operator A k i 
in the polar coordinates, we see that it commutes with the operator d ipi . for any l G C. 
Then, in view of the relation Jf £S = 2d Vk V k ff s established in Proposition 3.3 (i), it suffices 
to find a solution <p k j of equation 

A kj cf kj = 2V™ S + C(T k , Tj), (3.39) 

for some constant C(T k ,Tj) (possibly depending on the temperatures 7},7fc), such that 
the function r) kj := d Pk (p k j satisfies the conditions (a)-(c) from Lemma 3.2. By the direct 
computation we get that the function 

4>kj —\uj u k \ — (Tj T T k )\uj u k | 

satisfies (3.39). Using that d, Pk ( | u 3 — u k | 2 ) = —2 (iu k ) ■ Uj, we find 

3 

Vkj = d Vk 4> kj = ~(iu k ) • Uj(\uj - u k \ 2 + 4(7) + T k )). (3.40) 

It is not difficult to see that the function rj k j satisfies the required conditions. Note also 
that 

J™ = 2 S «V™ = -3 (iu k ) ■ - u k \ 2 . (3.41) 

Substituting (3.40) and (3.41) to the first representation for the conductivity k from 
(3.10), passing to the coordinates r-j := \uj\, r k := \u k \, ip — <fj — <Pk, and using that 

\ u j — u k\ 2 = f k + r j — 2 rkTj cosip and (iu k ) ■ Uj = r k rj sin ip, 

we get the desired result by the direct computation. □ 


3.3 Space-time correlations of the energy flow 


In this section we study stationary space-time correlations of the energy flow under the 
limits £, A —y 0 and investigate their connection with the conductivity k. Since till the 
end of the section we will not talk about the uniformity of convergences with respect to 
the choice of the set C, further on constants are permitted to depend on it. 


Theorem 3.12. Let be a Markov semigroup associated with eq. (2.25) and p £,x be 
its unique stationary measure. Take any k,j , m, l G C satisfying \k — j\ — \m — l\ — 1 and 
denote 

POO 

W’ A := / (//’\ J k 3 Wr X J m i) - (k £ ’\ Jkj)(T £ ’\ Jmi) dr. 

Jo 


Then 


lim lim y £ ’ x = 
A—^0 £—^0 


7j7fc/<7),7fc) 

—TjT k K(7j,T k ) 

0 


if m = k, l = j, 
if m = jj = k, 


(3.42) 


Corollary 3.13. Assertion of Theorem 1.8 is satisfied. 


Proof of Corollary 3.13. Let us first note that for the constant temperature profile 
T n = T Vn G C we have (//’ A , J k f) = 0, for any k, j G C, \k — j\ — 1. Indeed, since in the 
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present case the measure /a e,x is just the Gibbs measure, /i e,A = Z~ 1 (J~)e~ H " dLeb, this 
follows from evenness of the Hamiltonian H u with respect to the variable p, and oddness 
of the energy flow Jkj with respect to it. 

Now let us rewrite the r.h.s. of (1.18) as 

T- ,y P?o ^ n o ( ' ’ ^n+^r <Jjj+ 1 ) dr+ / (l l 5 TT/j+ilV Jkk+i) dr^j. 

1 3=0 o 0<3,k<N-l,rfk J 0 

Due to Theorem 3.12, the limits of the first summand in the brackets equal to lVT 2 k(T), 
while the limits of the second one vanish. □ 

Proof of Theorem 3.12. The proof of the theorem is based on the following auxiliary 
lemma, which we will establish in the end of the section. Denote by *jl A and the Markov 
semigroups associated with the effective equation (2.37) and the Ornstcin-Uhlenbeck pro¬ 
cess (3.1) correspondingly. Let 

oo oo 

:= j <VpJ ml (u) - <//'\ J mI ) dr, TiUu) := J <V x T J%‘(u) - {/, J%) dr, 

0 0 

oo 

-=-J q 

0 

where p x is the unique stationary measure of eq. (2.37). Note that, due to (3.5), the 
function r] m i coincides with that from Theorem 3.4. Put also 

OO 

€?(“):= J - (ii’- x ,J ml )\dT, 

o 

and define fj^u) and i) m z(u) in a similar way. 

Lemma 3.14. (i) There exist constants C,n > 0, independent from e and X, such that 
for any u E we have 

I VmlHlvtlH, \VmlHl Vml(u) < C{1 + |ti|”). 

(ii) For any u E we have 

-4- rj^fu) as e^O and rj^u) -4 -r) ml (u) as A -)• 0, 

where the both convergences hold uniformly in u E K, if K C is a compact set. 

(Hi) For any 9 E [0, 2tt) we have r)£ ll (e l0 u) = rj^^u). 

Let us first show that the integral y e,x converges. Due to Fubini’s theorem, for this 
purpose it suffices to establish that 

OO 

(p e ’\ J UkjiW/Jml - <//’ A , Jml))\dr) < oo. (3.43) 

0 

20 Recall that (fjP, <J^f) = 0, due to Proposition 3.3. zz, where /z° is the unique stationary measure of the Ornstein- 
Uhlenbeck process (3.1). 
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Note that the l.h.s. of (3.43) equals to (/P’ A , \Jkj\Vmi)- Since the energy flow J k j and 
the function rj £ m A have at most a polynomial growth, inequality (3.43) follows from the 
estimate for the stationary measure p £,x provided by Lemma 2.1 .ii. Now let us pass to 


the limit e —> 0. Fubini’s theorem implies that 

= (3.44) 

Theorem 2.8 . % joined with Lemma 3.14 ,i,ii implies 

(n £ ’ X , JkjVml) (/A JkjVml) as £ 0- (3.45) 

In view of Proposition 2.12 and Lemma 3.14.ra, we have 

(h A ; ^kjVml) = (/A ('-fkjVml) It) = (/A = (A> AA Vml) ■ (3.46) 

Now we pass to the limit A —>■ 0. Proposition 3. l.ii joined with Lemma 3.14 .ii implies 

as A —» 0. (3.47) 

Combining (3.44)-(3.47), we see that it remains to establish the identity 

{ TjTi,K{T r Tk) if m = k,l = j, 

-TjT,,K{T r T t ,) if m = j,l = k, (3.48) 

0 if {mj}^{k,j}. 


The first equality of (3.48) immediately follows from the formula (3.10) while the second 
one from the relation rjjk = — rjkj. The third equality follows from the rotation invariance 
of the Gaussian measure /i° and can be obtained by the argument used in the proof of 
Proposition 3.7. □ 

Proof of Lemma 3.1 f. Let us start with several auxiliary propositions. 

Proposition 3.15. The effective equation (2.37) is exponentially mixing with the rate 
which is uniform in 0 < A < 1 (in the sense as in Theorem A.l.) 

Proof. It can be shown by an argument similar to that used in the proof of Theorem A.l, 
but simplified. Note only that in the proof of the reccurence one should apply the Ito 

formula to the function e Sr Y) Ij- 21 D 

tec 

Proposition 3.16. For any compact set K C the time-integral from the definition 
of the function converges uniformly in u G K and e, while that from the definition of 
the function rj^ converges uniformly in u G K and A. 

Proof. Recall that in Appendix A it is proven that eq. (2.25) is exponentially mixing 
with the rate which is uniform in e. Then it can be shown that for all u £ we have 


m X Jmi(u) - (p. £ ’\jmi)\ < C(l + \u\ n )e~ bT , (3.49) 


21 This is related to the fact that the sum of actions h ' s a- first integral of the resonant Hamiltonian H res . 

j'6C 


latter follows from the relation { £ Ij,H res } = - £ d v H rBS = - £ £ d Vj V™° = - £ d Vj V™- 

iac jec jeCk-.\j-k\=i \i-k\=i 


Proposition C.I.to. 


The 

and 
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where the constants C,n,b > 0 are independent from e (see e.g. a proof of Theorem 1 from 

OO 

[PV]). It follows that the integral f ty £ T ' x J m i{u) — (//’ A , Jmi) dr converges uniformly in £ 

o 

OO 

and u G K. To prove that the integral J *p x J(%j s (u) — (/i A , JfJ) dr converges uniformly 

o 

in A and u G K, we use Proposition 3.15 and argue similarly. □ 

Now we turn to the proof of the lemma. 

Item (i). The estimates |^ A |, r )(f X < C(l + |u| n ) follow from (3.49). The other estimates 
can be obtained similarly. 

Item (ii). Let us first show that r/^ A (u) —v 77 ^(w) as e —* 0 uniformly in u G K. Due 
to Proposition 3.16, it suffices to prove that for any T > 0 we have 

T T 

J W x J ml (u) - <//’ A , J rnl ) dr^ J ^ X r JJf(u) - (/i\ JZ S ) dr as e —>• 0, 

0 0 

uniformly in u G A". Since, due to Theorem 2.8. * joined with Proposition 2.12, we have 
(// ,A , Jmi) —y (/i A , ) as £ —* 0, it suffices to establish that 

T T 

J ^ x Jmi{u) dr —> y* fp A ^® s (ii) dr as e —> 0 uniformly in u G AT. (3.50) 
0 0 

T 

Let us rewrite the l.li.s. of (3.50) as J E J m i(u £,x (r)) dr, where u £,x {r ) is a solution of eq. 

0 

(2.25) satisfying n e,A (0) = u. Denote by a E)X the process u £,x written in the a-variables. 

Remark 3.17. Simple analysis of proofs of Lemma f.l and Theorem 2.f shows that 
convergences obtained there for the process a £,x hold uniformly in u G K, if the set C is 
fixed. 

Applying Lemma 4.1 to the function h(u,6) := J m i(e l0 u ), in view of Remark 3.17 we 

T 

obtain that the integral above is uniformly in u G AT close to J E JJi s {cl £,X (t)) dr. Due to 

0 

Theorem 2.4, the latter integral is uniformly in u G K close to 

T T 

/ EjsVM)<fr= f <p r A j,:r(u)d T , 

0 0 

where u x {r ) is a solution of the effective equation (2.37) satisfying rz A (0) = u. 

The convergence 77 ^ (w) ~Vmi{u) as A —» 0 can be established with help of Proposi¬ 
tion 3.1 in a similar way. 

Item (Hi). Let u x {r) be a solution of the effective equation (2.37) satisfying u A (0) = 
u. The rotation invariance of the effective equation (see Proposition 2.11) implies that 
e l9 u x {r ) is its weak solution. Since it satisfies e l6 u x { 0) = e™u, we have 


rfi u (e*u) = j E JZ s (e l “u x (T))-( f r x ,J( n r)dr= J E JJf) dr = r£Ju), 

0 0 

where we have used Proposition C.l.i □ 
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4 Proofs of the averaging theorems 

Here we prove Theorems 2.4, 2.8 and Proposition 4.9 which implies Proposition 2.7. 


4.1 Averaging lemma 

We start with the following averaging lemma which is the main tools in the proofs of 
Theorems 2.4 and 2.8. Let u £ (t ) be a solution of eq. (2.25), either satisfying T>{u £ {fS)) = 
V{uq) or a stationary one. Let a e (r) be the corresponding a-variables, given by (2.27). 
Fix a bounded set A C Z d and assume that the set C satisfies C 3 A. Take a function 
h(u,9) G £, 0 C (C |c| x [0, 27 t)) such that for each 9 the function h(-, 9) satisfies Supp h(-, 9) C 
A and has at most a polynomial growth at infinity which is uniform in 9. Then there exist 
non-decreasing functions c, C\ : M i —> M such that 

\h(ui,9)\ < c(max |i/ij|), |/i(ui, 9) - h(u 2 , 0)| < c x ( max \u kj \) V' \u Xj - u 2j \, (4.1) 

j€ A k=l,2,j£A ' L —* 

jeA 

for any U\,U 2 G C^, 9 G [0, 2tt), where the function c has at most a polynomial growth 
at infinity. Denote 

2t r 

( h)e(u ) := j h(u,9)d9. 

o 

Lemma 4.1. Take a function h as above. Then 


E max 

0 <t<T 


(, h(a e (s),£ x s) — ( h)o(a £ (s ))) ds 


—y 0 as £ —y 0, 


uniformly in C and with respect to the choice of the function h satisfying (f.l ) for fixed 

C,Ci. 


Proof. For purposes of the proof we first introduce some notations. For an event P and 
a random variable £ we denote 

E r £ := E (£%). 

By ae(r) we denote various functions of r G 1 such that ae(r) — > 0 as r — y oo uniformly 
in C and h, satisfying (4.1) for fixed c, c\. By aeoo(^) we denote functions ae(r) such that 
ae(r) = o{r~ m ) uniformly in C and h. for each m > 0. We write ae(r) = ae(r; a) to indicate 
that ae(r) depends on a parameter a. 

Denote by Pi (A) the neighbourhood of radius 1 of A in C, i.e. 

Pi (A) := { n G C I there exists k G A satisfying \n — k\ < 1}. 


Fix R > 0. Set 

Q R = { max maxJ4 (r) | > R}. 

fce(7i(A) 0<r<T 

Due to the estimate (i) of Lemma 2.1, we have 


(4.2) 


P (Q h ) < teooiR). 


(4.3) 
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In view of the polynomial growth of the function h, 


E o max 

R 0<t<T 


h(a £ {s),e 1 s)ds 


— 3Soo(-^)> 


and the function ( h)g(a £ (s )) satisfies a similar relation. Then, it is sufficient to show that 
for any R> 0 


U £ := Eq max 

H 0 <t<T 


h(a £ (s),E 1 s) — (h) e (a £ (s)) ds 


—> 0 as £ —> 0 uniformly in C, h. 


For this purpose we consider a partition of the interval [0, T] to subintervals of the length 
6 = y/e by the points Ti = 15, 0 < l < L, where L = [T/5\. Denote 


n+i 


Vi — h(a £ (s),e 1 s) - (h) e (a £ (s)) ds. 


Then 


We have 


L—l 


U E <En R J2h\ + Mc{R). 


(4.4) 


/=o 



n+i 

p 


n+i 

p 

1 < 

/ h(a £ (s),e~ 1 s) — h(a £ (ri),£~ 1 s) ds 

+ 

/ h{a £ (r{),£~ l s) - (h) e {a £ (Ti)) ds 


U 

n 


n 


n +1 


(,h) e (a £ {ri )) - ( h)g(a £ (s )) ds 


n 

1 "\;3 


=-.y{ + yt + yf- 


Terms yl,yf. Since the function h is locally Lipschitz, 

{y} + yf) < 25d(R) y^Eq fl max |d(s) - aUn) \ < R), 

U Ti - s - Ti+i 

where the last inequality follows from eq. (2.29). 

Term y 2 . Changing the variable of integration s := e _1 s, we obtain 


(4.5) 


E n R y 2 — sEf 


£ 1 U +1 


h(a e (Ti),s ) - ( h)e{a £ (Ti))ds 


< 2 ec{R), 


(4.6) 


due to the definition of the function (h)g. 

In view of the identity e = 5 2 , estimates (4.5) and (4.6) imply 


E n *M <5( as (d- 1 -,R) + 26c(R)). 
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Noting that L <T5 1 , from (4.4) we get 

U £ < T(ae(d _1 ; R) + 2 Sc(R)) + 2 5c(R) —» 0 as £ —» 0 uniformly in C, h. 

□ 

Let now u £ {t ) be a stationary solution of eq. (2.25). 

Corollary 4.2. Let a function h G £/ oc (C^) satisfies Supp/i C A and has at most a 
polynomial growth at infinity. Then for any t > 0 we have 

E/i(« £ (r))-E/i(o £ (r))G0 as e —* 0 uniformly in C. 

Proof. There exist functions c, C\ : Mg R, where c has at most a polynomial growth 
at infinity, such that 

|h(«i)| < c(max|'Uij|), \h(uf) - h(u 2 )\ < ci( max \u kj \) V' \uij - u 2 j\, (4.7) 

j sA fc=l,2,jsA LJ 

j£ A 

for any Ui,U 2 G Cl c L Consider the function h(u,9 ) := h(e z9 u). Clearly, it satisfies (4.1) 
with c, Ci from (4.7). Applying Lemma 4.1, we get 

T 

j E/r(n e (r)) — E (h)i?(a e (r)) dr —?• 0 as £ —)■ 0, (4.8) 

o 

uniformly in C and with respect to the choice of the function h satisfying (4.7). In view 
of Proposition C.l.z, we have (h)ii(a £ (r)) = (h)/j(M e (r)). Then stationarity of the process 
u £ implies that the function under the integral in (4.8) is independent from time, so that 
for any r > 0 we have 

E h{u £ (r)) — E ( h) R (a £ (r )) g 0 as £ —)■ 0 uniformly in C and h. (4.9) 
Let h £,T {u ) := h(e~ l£ lr u). Since the functions h £,r satisfy (4.7) for any £ and r, we have 
E h £ ' T {u £ {r)) — E (h £,T ) r{o £ (t)) g 0 as £ —> 0 uniformly in C. 

Since h £,T (u £ (r )) = h(a £ (r )) and ( h £,T )n = (h)/j, we get 

E/i(a £ (r)) — E (h) R (a £ (r)) -GO as £ —> 0 uniformly in C. 

Jointly with (4.9) this implies the desired convergence. □ 

4.2 Proof of Theorem 2.4 

To prove the theorem it suffices to fold weak solutions a £ (r) and u(t) of eq. (2.29) and 
the effective equation (2.31) correspondingly, defined on the same probability space and 
satisfying X>(a £ (0)) = V{u{ 0)) = T>(u 0 ), such that for every j G C we have 

E max \w £ (t )\—)• 0 as £ —y 0 uniformly in C, (4-10) 

where w £ {r) := a £ (r) —u(t). Let us start with constructing suitable Brownian motions /? 
and (3. Denote 

T 

0‘ = (/3JW, /3J(r) := V2 J 

0 
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Proposition 4.3. We haveV{(3 £ {-)) ->■ X?(/3(-)) as e ^ 0 m C{[0,T],C ]Cl )■ 

Proof of the proposition is given after the end of the proof of the theorem. Due to 
Proposition 4.3 and Skorokhod theorem, we can find random processes f3 £ (r) and /3(r), 
defined on the same probability space and satisfying X?(/3 e (-)) = V(/3 £ (-)), V(/3(-)) = 
V(f3(-)), such that 

(3 £ —» (3 as e —» 0 in C([0, T], C' c ') a.s. (4.11) 


Take solutions a £ {r ) and u{r) of eq. (2.29) and the effective equation (2.31) on the 
probability space above (extended if needed), corresponding to the Brownian motions 


(3(t) := 


V2. 


e* £ ls df3 £ and (3, and having the same initial conditions, distributed as u$. 


Let us show that the process w £ = a £ — u satisfies (4.10). Due to eq. (2.29) and (2.31), 


T 


T 


K( r )l ^ 




S a £ (s )) — 1Zj(a £ (s)) ds 


+ 


lZj(a £ (s)) — 7 Zj(u(s)) ds 


o o 

+ v^l 0j(T) - ft(r)| =: y[j{T) + y^T) + yj,(r). (4.12) 

First we will estimate terms 3^,3^ and 3 7 | 7 - separately. 

Terms 3 7 f,- and yL. Using that E max \(3 £ (t)\ 2 V |/3Ur )| 2 < C, it is not difficult to 

J J 0 <t<T j 

show that (4.11) implies 



as e —» 0 . 


Since each of the vectors /3 £ and (3 has independent components, the processes /3J, (3j can 
be chosen independently from the choice of the set C, so that the convergence above holds 
uniformly in C. 

Applying Lemma 4.1 to the function h(u,9) = e~ ld Pj(e l9 u), we have 
E max yh(r) —> 0 as e —> 0 uniformly in C. 

0 <r<T 3 

Denote 5 £ = ( Sj)jec, where S £ E max^iy^iT) + y^ij)), so that 

S £ 0 as £ —> 0 uniformly in C. (4-13) 

Term yf 2j . Since, by assumption HV, the second derivatives of the potential V have 
at most a quadratic growth, the same holds for the first derivatives of the functions 7 Zj. 
Then for each j E C and any r > 0 we have 

T 

ytj <c f (1 + (i a fei 2 + KD) \ w i\ ds ■ ( 4 - 14 ) 

0 k:\k-j\<l k:\k-j\<l 

Now we go back to the proof of (4.10). Define for R > 0 a stopping time 22 
t r = inf{r > 0 : 3m G C satisfying |a^(r)| 2 V |M m (r)| 2 > R(\m\ 1 ^ 2 + 1)}. ( 4 -15) 

22 We can not define a stopping time in a standard way tr = inf{r > 0 : max|a^(r)| 2 V |w m (r)| 2 > R} since P (tr < 

m€C 

T)^0asJ?^-oo non-uniformly in C. To overcome this difficulty we introduce a stopping time tr which admits a slow 
growth of the absolute values |a^J, |u m | with respect to \m\. 
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For L e N denote \w £ \l ^ e Since for any s < r R and k e C satisfying 

!tl<i 

| k | < L + 1 we have 

K(s)| 2 V Ms)| 2 < R((L + 1) 1/2 + 1) < C(R)VL, (4.16) 

estimates (4.12) and (4.14) imply that for 0 < r < T we have 

tAt r 

Cj2 e ~ blE f (i+ (l a fc| 2 + K| 2 )) \ w k\ ds 

\j\<L o k:\k-j\<l k:\k-j\<l 

tAt r 

Ci(r)Vl e / E E e ^\w £ k \ds 

o lil<£fc:|fc-il<i 

tAt r 

C 2 (R)VLE f (\w £ \L + e~ L ^ |ta| ids 

0 |/c|=L+l 

T 

C 2 {R)Vl I (e\w £ (s Ar R )\ L +e~ L ^ \ w k (sAr R )\)ds. 

0 |fe|=£+l 

(4.17) 

Due to (4.16), for k G C satisfying \k\ = L + l we have \w k (s At r )\ < |a|(sAri?)| + |wfc(sA 
tr)\ < C(R)L 1 ^ 4 . Since the cardinality of the set {k : \k\ = L + 1} is bounded by CL d ~ 1 , 
we get 

E K(sAr*)|<C(fl)L“- 3/4 . 

|fc|=i+l 

By the Gronwall lemma, relation (4.17) joined with the inequality above implies 
E | w £ (t A t r )\ l < (|<f | L + C(R)e- L L d -*/ A )e c{R)VZr . 

In view of estimate (4.13), we have \S £ \l —> 0 as e —> 0 uniformly in C, for each L e N. 
Letting L —> oo and e —> 0 in such a way that \5 £ \ R e c ( R ' , '^' T —> 0, we obtain that 
E | w £ (t A t r )\l -A 0 uniformly in C. In particular, for each k G C and any 0 < r < T we 
have 

E | w £ k {r A t r ) | —y 0 as e —> 0 uniformly in C. (4-18) 

Inequality (4.12) jointly with estimate (4.14) implies 

T 

E \w £ (t A tr)| < 5 £ + C(R)(y/\j\ + 1) / E ^ \w £ k (s A t r )\ ds. (4.19) 

0 k:\k-j\<l 

Since the function under the time integral in (4.19) is bounded by a constant C(R,j), in 
view of (4.18) the dominated convergence theorem implies that this integral tends to zero 
as £ —> 0, uniformly in C. Then in view of (4.13) for every j E C we have 

E rriax^ \w £ (r A r R )\ —> 0 as £ —> 0 uniformly in C. (4.20) 

To recover convergence (4.10) from (4.20) we will need the following proposition. 


E | w £ (t A t r )\ l < \S £ \ l + 

<m L + 

< |5 £ |l + 

< \S £ \l + 
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Proposition 4.4. (i) P (tr < T) — y 0 as R — > oo uniformly in C. 

(ii) For all j G C we have E max e"l u A r )l < C , where a,C > 0 are the constants from 

[ 7 0<r<T ~ 

Lemma 2.1. 


Using the Cauchy-Schwarz inequality, we find 


E max \w £ At) \ — El, 


0 <t<T 


(• tr>t) max 

v ' 0 <t<T 


I Wj(T A Tr)| + EI (rfl < T ) nm |w £ (r)| 

1 / 2 , 


< E max |ud(r A 7r)| + (P (tr < T)) “(E max |w; £ (t)| 2 ) 

0<T<T J x / \ 1 4 x 7 1 / 


| 2 \ 1/2 


0<t<T 


Proposition 4.4. ii joined with the estimate of Lemma 2.1. i implies that E max \w e Ar)\' 2 <C 


0 <t<T 


Then, in view of (4.20) and Proposition 4.4. i, we arrive at (4.10). 

Proof of Proposition f.f. (i) We have 

P (t r < T) = P (3 j E C : max |a £ (r )| 2 v |u i (r )| 2 > i?(|j | 1/2 + 1)) 

= P (3 j eC: max |a £ (r A t r ) | 2 V \ Uj {r A t r ) | 2 > R(\j\ 1/2 + 1)) 

< E P ( A r ^)| 2 ^ ^(l^'l 172 + !)) 


□ 




p ( 0 ^ t <t \ Uj ( T A Ti? )l 2 - R (\i \ 1/2 + x ))- 

tec 


In view of the estimate E max e "l a j( TAT «)l 2 < E max e a ^ T ^ 2 < C provided by Lemma 2.1, 

0<r<T 0 <t<T 

convergence (4.20) joined with the Fatou lemma implies 


E max e “K ,(rAr fl)l 2 < c, je C. 

0<t<T ~ 


(4.21) 


Then, using the Chebyshev inequality, we obtain 


P (tr < T) < C e~ aR{ljll/2+1) <CJ2 e" QH(li|1/2+1) ^0 as R -> oo uniformly in C. 

j&C j£Z d 

(ii) The desired estimate follows from item (z) joined with (4.21) and the Fatou lemma. □ 
Proof of Proposition f.3. Since components of each of the vectors (3 and (3 e are inde¬ 
pendent, it suffices to prove that D(/3j(-)) — 1 T>(/3j(-)) as e — * 0 in U([0,T],C) for each 
j e C. In a standard way it can be shown that the set of measures (X>(/3j(-)), 0 < £ < 1 } is 
tight. Let Q be its limiting point as £*, —> 0. Take a process 0j(r) such that V(0j(-)) = Q. 
We need to show that / 3j is a standard complex Brownian motion. Since the processes f3) 

are uniformly in £ square-integrable martingales, the process (3j also is. Then it suffices 
to establish that 

\Pi]t = [/3 2 ]t = t and [/3i, j3 2 ] T = 0, 

where 0i = Re (3j, /? 2 = Im f3j, [-] r denotes the quadratic variation and [•, -] T stands for the 
cross-variation. We will only show that [0i\ T = t, proofs of the other assertions are similar. 
Due to the definition of the quadratic variation, for this purpose it suffices to prove that 
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the process 7 (r) := (/3i(t)) 2 — r is a martingale. Denote Y k ( T ) : = 
f5[ k := Re/3 £k . Clearly, we have 

D( 7 £fc (-))P( 7 (-)) as k —> 00 . 


On the other hand, since [fd[ k ] T = 2 f cos 2 (e k 1 s) ds, 


7 £k (r) = M £k (r ) + 2 / cos^(e fc 1 s) ds — r, 


(/3^(t )) 2 — r, where 

(4.22) 

(4.23) 


0 

where M £k (r) is a square-integrable martingale. Since max 


T 


2 / cos 2 (e fc 1 s) ds — r | —>■ 0 
0 


as /c —?• 00 , (4.23) jointly with (4.22) implies that 

D(M £ ‘(-))-D( 7 (-)) as fc^oo, 

so that 7 (r) is a square-integrable martingale as well. □ 


4.3 Proof of Theorem 2.8. 

Item (i). Let u £ {t ) and u{t ) be stationary solutions of eq. (2.25) and the effective equation 
(2.31) correspondingly, V{u £ {t )) = // and D(n(r)) = /i. Take a function h G £&(Cl c l). It 
suffices to show that 

E/i(!i £ (r))GE%(r)) as £ -> 0. (4.24) 

Let us pass to the a-variables a £ (r), corresponding to the process u £ (r). Consider the 
space C*([0, 00 ), C^) provided with the topology of uniform convergence on finite time 
intervals, which is given by the metrics 

OO 

p(v,w) = 'S~' 2~ k sup \v(r) — w(t)\ A 1. (4-25) 

0<t<K 

Proposition 4.5. The set of laws V(a £ (-)), 0 < £ < £ 0 , 23 is tight in C([0, 00 ), Cl c l). 

Proof. It suffices to establish that the set of measures T>(a £ (-)) is tight in C([0, T], C^) 
for any T > 0. The latter follows in a standard way from the estimates of Lemma 2.1, 
equation (2.29) and Arzela-Ascoli theorem. □ 

Consider some limiting point 

V(a £k (-))-+Q as e k 0. (4.26) 

It turns out that the measure Q is a law of a weak solution of the effective equation: 

Proposition 4.6. The measure Q coincides with the law T>(u°(-)) in the space C([ 0, oo),Cl c l), 
where u°(r) is a weak solution of the effective equation (2.31). 

23 Recall that £q is defined in Lemma 2.1. 
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Proposition 4.6 can be established by argument similar to that used in the proof of 
Theorem 2.4. 24 We do not prove it here since below we will establish an analogous result 
in a more complicated, infinite-dimensional setting (see Proposition 4.8). Proposition 4.6 
provides that for any r > 0 we have E h(a £k {r)) —> E h(u°(r)) as e k —> 0. Jointly with 
Corollary 4.2 this implies 

Eh(u £k (T)) ->• E/i(u°(r)) as e k -> 0. (4.27) 

Since the process u £ is stationary, (4.27) implies that the process u°(t) also is, so that 
V(u°(t)) = V(u(t)) = n . Thus, we get (4.24). 

Item (ii). Assume that condition HC a0 is satished. Recall that the uniformity in C 
of the weak convergences of measures through all the text is understood in the sense of 
finite-dimensional projections. That is, we need to prove that for any bounded set AcZ d 
and a function h G £&((d c l) satisfying Supp h C A, convergence (4.24) holds uniformly in 
C satisfying C 3 A. For this purpose it suffices to show that for any sequence of bounded 
sets (Cn)neN, C n C satisfying C n D A for each n G N, there exists a subsequence (C n )neN 
such that the desired convergence holds uniformly in C G (C n )neN- 

The proof follows a scheme below. As (C n ) ne ^ we take an arbitrary subsequence which 
has a "limit" as n —> oo, say Coo- Ass. HCoo provides that the corresponding Coo-effective 
equation has a unique stationary measure in A4 Co °. Then, arguing as in item (i), we 
establish convergence (4.32), where u £,Cn denotes a stationary solution of eq. (2.25) with 
C = C n , while u°° denotes that of the Coo-effective equation. This provides control of 
the expectations E h(u £,Cn (r)) for large n and small e, sufficient to prove the desired 
uniformity of convergence. Now let us present the rigorous proof. Denote by B N a ball 
in (Z d , | • |) of the radius N centred at zero, where | ■ | denotes the 1 1 -norm in 7L d . 

Proposition 4.7. Let (C n ) ne pj be a sequence of bounded sets in 7L d . Then there exists a 
subsequence (C n )neN and a set Coo C Z d such that C n —> Coo as n —> oo in the sense that 
for any N G N and all n > N we have C n D B^ = Coo H Bjy. 

Proof. Since the number of sites in Bi is finite, there exists a subsequence (C nni ) niS M 
such that for every j G B i we have either j G C n or j ({ C n at the same time 
for all n\. Then the set C n fl B\ is independent from n\. Similarly, we hnd a sub¬ 
sequence (C nni )n 2 eN such that the set C„ ni fl B 2 is independent from n 2 - We con¬ 
tinue the procedure and apply the diagonal process to choose a subsequence (C n ) n£ N = 
C ni ,C n , C n ,.... Then for any N G N the sets C n fl B N coincide for all n > N. 

Let Coo '■= U neN (C n fl B n ). Clearly, the pair (C n ) nS N and C a0 satishes the assertion of the 
proposition. □ 

Choose a subsequence (C n )neN of the sequence (C n )„ S N as in Proposition 4.7. Let 
u £ ' Cn (r) be a stationary solution of eq. (2.25) with C = C n , and u Cn {r) be that of the 
effective equation (2.31). We need to prove that convergence (4.24) holds uniformly in 
C G (C n ) nS N, i.e. to show that for any 5 > 0 there exists e > 0 such that for any e < e 
and all n G N we have 

|E h{u e ^ n {r)) — E h(u^ n (r))\ < 6. (4.28) 

24 Note that this argument does not lead to the uniformity in C of convergence (4.26). Indeed, in the present case the 
distributions of initial conditions T>(a £k ( 0)) are different for different so we should add in (4.17) the term E \w £k (0)|l- 
We do not know if it convergences to zero as —»■ 0 uniformly in C. 


40 



Denote u £,n = (u £ f n ) , where 


u j ■= 


u 


£,Cn 
j ’ 
0, 


if jeCooDCn, 
if j e Coo \ C n . 


(4.29) 


Let a £,n = (0'j >n )j£c oo ^ ie corresponding a-variables. Consider the space C([0, oo), Cl c °°l) 
provided with the topology of uniform convergence on finite time intervals which is given 
by the metrics (4.25), where the distance | v(t) — w(t) | is replaced by 

P c \Coo\Ht),w(t)) = 2~ bl \Vj(r) — Wj (r)| A 1. 

j £Coo 

Using the uniformity in C of estimates from Lemma 2.1 we get that the set of laws 
{V(a £,n (-)), 0 < £ < 1, n 6 N} is tight in the space C([0, oo), C^°°l). Let Q°° be its 
limiting point as £*, —* 0, rik —> oo. 

Proposition 4.8. In the space C([0, oo), Cl c °°l) the measure Q°° coincides with the law 

U(u°°(-)), where u°°(t ) is a weak solution of the Coo-effective equation. For any j e Coo, 

0 < A < 1 and r > 0 it satisfies 

E max e a|u r (s)| < C. (4.30) 

se[r,r+l] 

Before presenting the proof of the proposition, let us finish the proof of the theorem. 
Since A C C n for all n 6 N, we have A C C^. Then, due to the inclusion Supp h C 
A, the function h can be considered as h : Cl c °°l i —> M. Proposition 4.8 implies that 
E h(a £k,nk (r)) —> E/i(n°°(r)) as £k —)■ 0 ,nk oo. Since, due to Corollary 4.2, we have 
E/j(r‘ A (r)) — E h(u £k ’ nk (r)) —> 0 as £k —> 0, —> oo, we get 

E/i(w £fc,nfc (r)) —>• E h(u°°(r)) as £k —> 0, nk —> oo, (4-31) 

so that u°° is a stationary process. Due to (4.30), the stationary distribution V(u°°(t)) 
belongs to A4 C °°. Then, in view of ass. HCoo, it is defined uniquely, so that the limit in 
(4.31) does not depend on the sequence ( £k,nk)keN and holds as e — > 0, n — > oo. Since 
Supp h C A C Coo n C n for all n E N, we have h{u £,n {r )) = h(n e ’ Cri (r)), so that 

E h{u £ ^ n {r)) —> E /i(tt°°(r)) as £ —» 0, n —> oo. (4.32) 

Due to item (i) of the theorem, for each n 6 N we also have 

E/i(w £,c ~ n (r)) —> E/i(u Cn (r)) as £ —> 0. (4.33) 

Jointly with (4.32) this implies that E h(u Cn {r)) —> Eh(n°°(r)) when n —> oo. Conse¬ 
quently, for any 8 > 0 there exist N E N and £\ > 0, such that for every n > N and 
0 < £ < eq, we have 

|Eh(n £ ’^(r)) — E/ i (n 0O (r))|, \Eh(u dn (r)) - E/i(tx°°(r))| < 8/2. 

Then for n and e as above, we get (4.28). In view of (4.33), we can choose e 2 > 0 such 
that (4.28) is also satisfied for every 0 < e < £2 and n < N . Then it holds for all n E N 
and £ < £1 A e 2 - Cl 
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Proof of Proposition f.8. Before starting the proof let us make the following remark. 
Since Supp P 3 , Supp Rj C {k G Tf : \k — j\ < 1} fl {k G C}, the functions 25 Pj and 7 Zj 
depend on the choice of the set C. To indicate this, we denote those corresponding to the 
set C nk as P^ IZj and those corresponding to the set as Pf°, 7l°°■ Now we start the 
proof. For shortness we write a k := a Ek ’ nk . Consider the process b k = (b k )j e ^, where 

T 

b){r) = a k (r) - aj(0) - J H™(a k (s)) ds, j G C*. 

o 

Consider also a process u°°(t) satisfying V(u°°(-)) = Q°°, and put b = ( bj )j e c 00 , where 

T 

bj (t ) = uf{r) - u?(0) - f K?(u°°(s)) ds , j G 4, 

o 

The convergence V(a k (-)) — 1 Q°° as k —> oo implies 

V(b k (-)) V(b{-)) as k^ oo on C([0, oo), C |doc| ). (4.34) 

On the other hand, according to (4.29) and eq. (2.29), for j G C ^ 0 C nk we have 

6j(r) = e}(T) + ^(r), (4-35) 

* „ „ _ T . _j 

where (3 k = {(3 k ) J - g c oo , /^(t) := ^/27j / (/ ®j)j & c oa are standard real indepen- 

oo o 00 

dent Brownian motions, and 

T 

0 J(t)= | (e- i£kls P k {e i£kls a ek ’^ k {s))-n°°{a k {s)))ds, 
o 

where a £k,Cn k denotes the a-variables corresponding to the process u £k,Cn k . Since C nk —> 4 
as k —> oo, for sufficiently large k (depending on j ) we have P k = Pf° and 7 Z™(a k (s)) = 

1Z c f(a £k,Cn k (s)). Then, applying Lemma 4.1 to the function h(u,6 ) = e~ l6 P°°(e ld u) , we 
get 

E max |0((r)| —> 0 as k — > oo for every T > 0. 

Moreover, Proposition 4.3 implies that V0 k (-)) —^ U(/3(-)) as k —» oo in C([0, oo), C^°°l), 
where {3 = (Pj)j e c 00 i Pj = \/7jPj, and (Pj)je.c x are standard complex independent 
Brownian motions. Then from (4.35) we obtain 

V(b k (-)) ^V0(-)) as k^ oo on C([0, oo), C 1 ^ 1 ). 

Thus, due to (4.34), we have T>(b(-)) = T>0(-)), so that u°°(r) is a weak solution of the 
Coo-effective equation. Estimate (4.30) follows from the Fatou lemma and Lemma 2.1. □ 

25 We recall that the functions Pj and IZj are defined in (2.26) and (2.32) correspondingly 
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4.4 The Coo-effective equation 

In this section we prove that under conditions of Proposition 2.7 the Coo-effective equation 
defines a Markov process and is mixing, in a suitable space. In particular, this will imply 
Proposition 2.7. Denote 

IMIo := ^l lkl \ u k\\ 

k&C 

where u G and C = C or C = Coo (the choice will be clear from the context). Fix the 
constant 1/2 < 7 < 1 as in Lemma 2.3. Consider the space 

X := {u G C |Coc| : ||w || 0 < 00}, 

and denote by V(X) a space of Borel probability measures on X. 

Proposition 4.9. Assume that conditions of Proposition 2.1 are fulfilled. Then for any 
unbounded set Coo Q ^ d the corresponding Coo-effective equation defines a Markov process 
in X and is mixing. That is, it has a unique stationary measure /x°° in the class of 
measures V{X), and for any its solution u{r) satisfying P (u(t) G X for all r > 0) = 1 
we have the convergence V(u(t )) — 1 pp° as r —* 00. 

I 11 order to prove Proposition 4.9 we will need the following 

Proposition 4.10. Let u(t) be a solution of the effective equation (2.31) satisfying 
E e <5olK°)ll 0 < C for some 5 0 > 0. Then there exists 8 > 0 such that for any 0 < A < 1 
and t > 0 we have 

E max e 5||u(s)ll ° < ^( 7 , 5, C). 

se[r,r+l] 

Proof. Note that the norm || • || 0 restricted to Rl c l coincides with the norm from (2.3) 
with j = 0. Let u = p+iq. Since ||w||q = ||p||o+IMIo> ^ ie u PP er bound from Proposition 2.2 
implies that Ee {olf °^^ 0 ^ < C for some 80 > 0, where U°(u ) := U°(p(u), q(u)). 26 Consider 
a solution u e (r ) of eq. (2.25) satisfying D{u £ { 0)) = T>(u( 0)). Due to Lemma 2.3, there 
exists 8 > 0 such that 

E max e ^° ( “ e(s)) < C u 

sG[t,t+ 1 ] 

for any 0 < c < £ (j . 0 < A < 1 and r > 0. Due to the lower bound from Proposition 2.2, 
we have U°(u ) > ||w||o/4, so that 

E max e 5|l “ £(s)l| o < C u 

se[r,r+l] 

where 8 := <5/4. Then Theorem 1.1 joined with the Fatou lemma implies the desired 
estimate. □ 

Proof of Proposition f.9. The proof is divided into three steps. In the first one we show 
that if A is sufficiently small, then the Coo-effective equation is contracting in the space 
X. In the second step we prove that the Coo-effective equation is well-posed and defines 
a Markov process in X. In the third one we show that the Coo-effective equation has a 
stationary measure from V(X), prove that the latter is unique, and obtain the desired 
convergence. 

26 See (2.2) for the definition of U°. 
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Step 1. Let w 1 (r) and u 2 {r) be two solutions of the Coo-effective equation having de¬ 
terministic initial conditions uI,Uq G X and corresponding to the same Brownian motion 
{/3j)jeC x - Let w := u 1 — u 2 . Since the second partial derivatives of the potential Vff s are 
bounded, for all j G Coo we have 


d | Wj 

dr 2 


<A C ^2 \ w k\\Wj 

k:\k-j\<l 



(4.36) 


see (2.37). Multiplying the both sides of (4.36) by 7 ^ and summing over j G Coo, we find 

-^\\w\\l < (XC - l)\\w\\ 2 0 , a.s. 

Applying the Gronwall inequality, we get 

IH'Ollo < lk(0)|loe (AC_1)r , a.s. 

Assume A < 1/C. Then 

\\w{t )\\ 0 —> 0 as r —» 00 , a.s. 


(4.37) 

(4.38) 


Step 2. Let us show first that the Coo-effective equation admits a weak solution in 
the space X. Let u/f G X be a deterministic initial data. Take a sequence of sets 
C n Coo H {j G Z rf : \j\ < n}. Let u Cn (r) be a solution of the effective equation (2.37) 
with C = C n and initial conditions Uq" G satisfying u qJ = for every j G C n . We 

define the process u n = (Uj)jeCooj where u™ = if j G C n , and u™ = 0, if j G Coo \ C n . 
Proposition 4.10 joined with the Arzela-Ascoli theorem implies 

Proposition 4.11. The set of measures {T>(u n (-)), n G N} is tight in C([0, 00 ), Cl c °°l). 
Let 

V(u nk {-)) Q°° as k —y 00 . (4.39) 

Take a process u°°(r) satisfying = Q°°. It is possible to show that u°°(r) 

is a weak solution of the Coo-effective equation. Obviously, it satisfies u°°(0) = u/f. 

Moreover, Proposition 4.10 joined with convergence (4.39) and Fatou’s lemma implies 

that E max e ,5 ll u00 0)llo < C, so P (u 0O (r) G X for all r > 0) = 1. 
se[r,r+l] 

Estimate (4.37) implies the pathwise uniqueness in X of solutions for the Coo-effective 
equation. Jointly with existence of a weak solution, by Yamada-Watanabe arguments (see 
[Yor, RSZ]), the pathwise uniqueness implies existence of a strong solution. Then, using 
arguments from Chapter 7 of [0ks], it can be shown that the set of solutions corresponding 
to all possible initial data from X form a Markov family. 

Step 3. Let & T (u,-) be the corresponding transition function. Convergence (4.38) 
implies that for any Mq,Mq G X we have 

||^V(no, •) - ^ t (u 2 0 , -)\\* Lip —^ 0 as r —>■ 00 , 

where || • \\* Lip denotes the dual-Lipschitz norm in V(X). It follows that, if the Coo-effective 
equation has a stationary measure p°° G V(X), then it is unique, and for any solution 
u(t) of the Coo-effective equation satisfying P (u(t) G X for all r > 0) = 1 we have 
T>(u(t)) —*■ /i°° (see Step 1 form the proof of Theorem 3.1.3 in [KS] ). To show that the 
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stationary measure exists, we employ an argument similar to that used in the beginning 
of Step 2 joined with the estimate of Proposition 2.5. ii. □ 

Now let us show that Proposition 4.9 implies Proposition 2.7. 

Proof of Proposition 2.7. Existence. Estimates of Proposition 2.5 imply that the 
stationary measure /i°° of the Coo-effective equation constructed in the proof of Proposi¬ 
tion 4.9 is regular. 

Uniqueness. Let p°° be a regular stationary measure of the Coo-effective equation. 
Then there exists a regular weak solution u°°(r) of the latter satisfying T>(u(t )) = p°°. 
Clearly, the regularity implies that P (u°°(t) £ X for all r > 0) = 1, so that p°° is 
a stationary measure of the Markov process in the space X, given by the Coo-effective 
equation. Proposition 4.9 provides that such measure is unique. □ 


A Appendix: Uniformity of mixing 


In this appendix we show that eq. (1.6) is exponentially mixing with uniform in £ rate. 
Note that for fixed £ the mixing property of eq. (1.6) is well understood, see e.g. [MSH], 
where it is proven for the case of smooth interaction potential V. It is convenient to work 
in the complex variables u = p + iq £ C^, so that instead of (1.6) we deal with eq. (2.25). 
The proof is based on the fact that the latter, written in the a-variables, has uniformly 
in £ bounded coefficients (see eq. (2.29)). 

We do not follow dependence of the rate of mixing on the choice of the set C, so in this 
section the constants C,C\,... are permitted to depend on C. 2 ‘ Moreover, for simplicity 
of notations we put A = 1. Let fJL; be a Markov semigroup associated with eq. (2.25). 

Theorem A.l. Equation (2.25) has a unique stationary measure p £ . There exist con¬ 
stants £i > 0 and C,b > 0 such that for all 0 < £ < £i and any Borel probability measure 
p satisfying ( p , \u\ 2 ) <00 we have 

\m*P - AIU < c( 1 + (p, M 2 ))e -br , (A.l) 

where || • || mr denotes the variational norm. 

Proof. In this proof by B R we denote a closed ball in of radius R, centred at zero. 
By 5P £ (u, •) we denote the transition function of the Markov process (2.25). 

A standard Bogolioubov-Krylov argument implies existence of a stationary measure p £ 
for eq. (2.25). Assume that for any U\,U 2 £ C^, 0 < £ < £1 and r > 0 we have 

\\^r(ui, ■) - ^ 2 , -)IU < C( 1 + Kl 2 + \u 2 \ 2 )e~ br . (A. 2 ) 


Using arguments from the introduction to Section 3, and of Sections 3.1.1,3.1.2 from [KS], 
it can be shown that (A. 2 ) implies the assertion of the theorem. Coupling argument from 
the introduction to Section 3 of [KS] implies that to prove (A. 2 ) it suffices to show that 
the Markov process given by eq. (2.25) satisfies the following two properties. 

Recurrence. Let u £ (t ) be a solution of eq. (2.25) satisfying u e (0) = u\, where u\ £ C^ c L 
Put 


t £ r \= inf{r > 0 : u £ (r ) £ B R }. 


27 We do not know if the rate of mixing is uniform with respect to the choice of the set C. If it is, it seems to be a very 
complicated problem. Note that we do not need this uniformity to get in Theorem 2.8 the uniformity of convergence in C. 
Indeed, in the proof of the latter we do not use the mixing property but we work directly with stationary solutions. 
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Then there exist R, S > 0 and £ 2 > 0 such that for any 0 < £ < £ 2 we have 

Ee k « < <7(1 + M 2 ). (A.3) 

Squeezing. For any R > 0 there exists constants 0 < d < 1 and £3 > 0 such that for 
any u\,u 2 € Br and 0 < e < £3 we have 

||^T(«1, •) - ^>{{u 2 , Ollt-or < 1 - 0- (A.4) 


We start with the proof of the reccurence. 

Proof of the recurrence. For simplicity of notations we will skip the upper index £. Let 

U(u) = H{u) + -^Pjqj, 

3 EC 

where u = p + iq, p = ( Pj)jec , q = ( qj)jec • Applying Ito’s formula to the function 
e ST U(u(r)) with 0 < <5 < 1, we obtain 28 

e 5 T U(u(r)) = U{ Ul ) + [ e 5 s ( 6 U + £~\U, H} - ^ Pj d Pj U + ^ TjcfaU) ds 

jec jeC 


^e Sa y/Wid^U dfa 

3 EC 


Using estimate (1.20), we get 

1 


Hu, h} - Y 1 Pk t <h } + I S v ( qm > 

Z j£C kec Z j&c Z Ifc—m|=l 


jec kec 

<5(IpI 2 -M 2 ) + Ci=(i + |9D. 


Moreover, for £ sufficiently small, 


\u\ 


V4 < u(u) < C 2 { 1 + M 2 ), - ^PjdpjU < (-1 + e)\p \ 2 + e\q\f 


3 ec 


Y. T - a lf = ^ 

3 EC j£C 

For K > 0 put £ := rjj A K. Then (A.5) joined with (A. 6 ) and (A.7) implies 
Ee*t/(u(£)) < U{uf) + E [ e Ss (A\u \ 2 + <7 4 ) ds, 


(A.5) 


(A. 6 ) 


(A.7) 


(A. 8 ) 


where 


A := [(hC' 2 +l/2-l+£)V(hC' 2 -l/2+C' 1 £+£)] = [(—l/2+hU 2 +£)V(—l/2+hC' 2 +£(U 1 +l))], 

and C 4 := C 2 + C 1 + C 3 . Choose £ 2 and 5 in such a way that for £ < £ 2 we have A(£, S) < 0, 
and fix R > -yj—C^/A{ e 2 , h). Let s < £. Then |w(s)| > R, so that A|m(s )| 2 + C 4 < 0, for 
any £ < £ 2 . Thus, (A. 8 ) implies 

E e s *U(u(£)) < U(u 1 ). 

28 Here it is simpler to look at eq. (1.6) than at (2.25). 
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Since, in view of the first estimate from (A.7), we have U(u(£)) > |«(£)| 2 /4 > R 2 / 4, we 
obtain 

Ee* < AU( Ul )/R 2 < 4C 2 (1 + \ui\ 2 )/R 2 . 

Letting K —y oo and applying Fatou’s lemma we get (A.3). □ 

Proof of the squeezing. Denote by PP a ' £ { y s, r; u, •) the transition function of the Markov 
process (2.29), corresponding to initial conditions u and initial time s. Let £Pf ,£ (u, •) : = 
^ a ’ £ (0, t; u, •). In view of the definition of the a-variables (2.27), we obviously have 

||^T(«l, •) - ^l(«2, ') \ \var = || («1, ') - 3*i e (u 2 , -) 11 «ar- (A.9) 

Let Q G B{ C^). First we will obtain uniform in e and u G B R estimates from below for 
the transition probability £Pf’ £ (u, Q). For this purpose we will use Girsanov’s theorem 
(see [0ks]). Let a £ (r) be a solution of eq. (2.29) satisfying a £ (0) = u, where u G B R , and 

= inf{r > 0 : |a £ (r)| > K}. 

In order to satisfy Novikov’s condition needed for application of Girsanov’s theorem (see 
below), instead of considering the process a £ (r) we will consider the process a £,A (r) which 
coincides with a £ (r) for r < and for t > satisfies 

af = ^2T je - i£ ~ lr j G C. 

Writing eq. (2.29) in the real coordinates and using that Pj are real valued functions, we 
see that Girsanov’s theorem implies that for all 0 < r < 1 the process a £,K {r ) satisfies 
the equation 

af = f)f , af(0)= Uj , j G C, (A.10) 

where (/ 3 £,K )j^c is a standard |C (-dimensional real Brownian motion with respect to the 
measure Q e,K given by 

dQ £ ’ K (co) = M £,K (cu)dP (cj) (A.11) 

with 

M £ ’ K = ex p(- / ^^(r,a £ A( T ))d^*| / |P £ (r, a £ ’*(r))| 2 dr), 

o o 

where by V £ = (Vj)jec we have denoted 

V’(t,u) := (27-)- 1/2 P J (e“"‘ T !1 ). 

1A P K 

Note that Novikov’s condition E exp J \V e (r, a £ ' K {r))\ 2 dr'j < oo, required for the 

o 

application of Girsanov’s theorem, is satisfied since for t < f e K we have |a £,A (r)| < K. 
Now let us estimate the transition probability &f’ £ {u,G) = P (a £ (l) G G) ■ We have 

P (a £ (l) G^)>P (a £ (l) G G, max |a £ (r)| < K ) (A.12) 

= P (a £ ’ K (l) G G, max |a £ (r)| < K) > P {a e ’ K { 1) G G) - a e(K), 
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where Lemma 2. l.i joined with (2.28) implies that a e(K) —> 0 as K —> oo uniformly in e 
sufficiently small. Due to (A.11), for any L > 0 we have 


P (a £ ’ K (l) eG)= I (M**)- 1 dQ e ’ K > e~ L Q £ ’ K (a £ ’ K (l) G G,M £ ’ K < e L ) 

{a e ’ K (1 )gG} 

> e~ L (Q e ’ K (a e ' K ( 1) EG)- e~ L ) , (A.13) 

where we have employed the exponential supermartingale inequality. In view of (A. 10), 
the random variable a e,K ( 1) has the Gaussian distribution with mean u, with respect to 
the measure Q e,K . Its covariance matrix written with respect to the real coordinates has 
/ D u D 12 \ 

the form D = diag I n { 2 r \22 ) , where 

V ■; J j^c 


Df = 27) / cos 2 (e 1 t) dr — 7)(1 + - sin( 2 e *)), 


1 

D f = 2Tj j sin 2 (e~ 1 r) dr = 7)(1 - | sin( 2 £ -1 )), 


D 1 j 2 = 2Tj I sin(e V) cos(e V) dr = : -^-(l — cos( 2 e 0). 


Since for e sufficiently small we have D > — diag 


Tj 0 

0 75 /^ 


, we hnd that 


Q 


e,K („e,K 


(1) eg) > Q £ ’ K (a £ ’ K (l) egn B x ) > C(u)Leb{g n Si), 


(A. 14) 


where Leb (g D Si) denotes the Lebesgue measure of the set g n B\. Obviously, 


C := min C(u) > 0 . 

uGBr 


(A. 15) 


Combining estimates (A. 12)-(A. 15), for any u e Br we get 

&>“’ £ (u,g) > e~ L (CLeb(gnB 1 ) -e~ L ) - ae(AT). (A. 16 ) 

Now we are able to estimate the variational norm (A.9). Without loss of generality we 
assume &l' e (ui,g) > (u 2 , g). Then, in view of (A.16), 

- &>r(u 2 ,g)\ = 1 - - &>r(u2,g) 

< 1 - e - i (C'Leb(Si) - 2e _i ) + 2ae(AT). 

Choosing hrst L and then K sufficiently large, we obtain 

\&>l’ £ {u l: g) - &?' £ (u 2 ,g)\ <1-0, where 0 := e~ L (CLeb{B 1 ) - 2 e" L ) - 2as(K) > 0 . 

It remains to note that the constant 0 is independent from e, u\, u 2 G Br and g G S(C^). 
Then, in view of (A.9), we arrive at (A.4). □ 
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B Appendix: Low temperature regime 

Let (p,q) = ( PjiQj)jec G M 2 I C L Consider the equation 

= Vh J t Pj = -<Li - £ - epj + J e c, 

k:\k-j\ = l 

where £, <5 1. It describes a system of weakly coupled oscillators, weakly interacting 

with thermostats of temperatures 87), j G C. Let us put 

Pi '■= V5pj, qj := VSqj, 

and assume that the interaction potential V is a homogeneous polynom of degree m > 3. 
Then 

j f Qi = Pi > j t Pi ~ ~<li ~ eS {m - 2)/2 J2 d <a V > 9k) - £ Pi + V^Tij/i > J e C - 

k-\k~i\=l 

Putting A := 8^ m ~ 2 ^ 2 , we arrive at equation ( 1 . 6 ), written in the fast time t. 

At the physical level of rigor the condition above for the form of the interaction poten¬ 
tial can be weaken. Since we study small amplitude solutions, the interaction potential 
V(qi,cfk) can be replaced by a leading order term of its Taylor series. Thus, it suffices to 
assume that the latter is a homogeneous polynom of degree m > 3. 

C Appendix: Resonant averaging 

Here we discuss some properties of the resonant averaging, given by formulas (2.33) and 
(2.34). Since we will use the derivatives with respect to the angles ip = (ipj)j G c, let us 
note that if a function / : t—R is C 1 -smooth, then / is continuously differentiable 

with respect to the angles (p. Indeed, this follows from the formula 

d Vj f = iUjd u .f - iujduJ. 

Proposition C.l. Let f : —> R be a continuous function. Then 

(i) For any £ G [0,27t) we have {f) R {e l ^u) = (f) R (u). 

(ii) Letf G A(C |C| ), Aoc(d c l) orC n (& c I), wheren G N. Then(f) R G £ 6 (C |C| ), A oc (d c ') 
or C n ( Cl c l) correspondingly. If partial derivatives of the function f have at a most poly¬ 
nomial growth at infinity, then those of (f) R also have at most the polynomial growth. 

(in) Let f G C' 1 (Cl c l). Then for every j G C we have {d tf , j f) R = d Vj {f) R . 

(iv) Let f G C 1 (Cl c l) be such that Supp / C {j , k}. Then d Vj (f) R = —d Vk (f) R . 

Proof. Items (i) and (ii) follow from formula (2.34). Item (Hi) follows from (2.33). To 
prove item (iv) it suffices to note that, in view of (2.33), the resonant averaging ( f) R (I , (p) 
depends on the angles pj, tp^ only through their difference <pj — (pk- □ 
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